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The  dynamic  instability  of  thin  shallow  spherical  shells 
under  uniformly  distributed  impulsive  loadings  has  been  discussed 
by  use  of  the  stability  theory  of  nonlinear  mechanics  together 
with  a  topological  method. 

The  axisymmetrical  buckling  of  the  shallow  shells  under 
three  types  of  dynamic  pressure,  namely,  an  impulse  function,  a 
step  function,  and  a  linear  function  of  time,  has  been  solved 
based  upon  a  first  approximation  and  Galerkin' s  method.   The 
results  of  the  first  two  cases,  within  the  range  of  the  validity 
of  the  theorem,  are  compared  respectively  with  those  obtained  by 
Humphreys  and  Bodner,  Budiansky  and  Roth,  and  Jiro  Suhara„ 

A  buckling  criterion  based  upon  the  relation  between  the 


dynamic  load- deformation  curve  and  its  counterpart  in  the  static 
case  has  been  established  through  the  comparison  of  the  results 
obtained  by  this  criterion  with  those  obtained  by  the  application 
of  the  stability  theorem  to  autonomous  systems. 

An  extension  of  the  same  techniques  to  the  solution  of 
problems  in  the  case  of  a  complete  spherical  shell  has  also  been 
made.   Based  upon  an  assumed  new  buckling  mechanism,  the  static 
buckling  pressure  for  such  a  shell  obtained  by  a  first  approx- 
imation and  axisymmetrical  deformation  theory,  agrees  well  with 
recent  experiments  conducted  individually  by  Krenzke  and 
Thompson   Transitional  instabilities  of  a  complete  spherical 
shell  under  uniformly  distributed  dynamic  pressures  in  the  form 
of  a  step  function  and  an  impulse  function  were  also  discussed. 
Results  are  presented  in  the  form  of  figures. 


INTRODUCTION 
1 .   A  Historical  Review  and  Recent  Advancement 

The  problem  of  dynamic  instability  of  a  thin  spherical 
shell  under  time  dependent  external  forces  is  inseparable  from  the 
problem  of  the  transverse  vibrations  of  such  a  structure.   Studies 
on  the  latter  problem,  as  in  the  works  of  Mathieu,  Lamb,  Lord 
Rayleigh  and  Love,  date  back  to  1882   As  may  be  found  either  in 
Love's  Theory  of  Elasticity  or  Rayleigh' s  The  Theory  of  Sound, 
these  classic  works  have  been  mainly  devoted  to  finding  the 
smallest  natural  frequency  or  the  "gravest  tone"  of  the  systemo 
The  method  used  by  Lord  Rayleigh  was  essentially  an  energy  method 
plus  an  assumed  displacement  pattern;  the  bending  and  membrane 
energies  were  considered  separately,  depending  on  whether  the 
middle  surface  of  the  shell  is  extensible  or  inextensibleo 

The  first  rational  study  of  the  transverse  vibration  of 
a  shallow  spherical  shell,  using  three  simultaneous  equations  of 
motion  for  the  three  displacement  components,  was  due  to  Feder- 
hofer  (44) „   Federhofer's  problem  was  discussed  later  by  Reissner 
using  a  different  approach  (45) ,   These  two  authors  were  interested 
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in  the  transverse  free  oscillations  of  a  shallow  shell.   In  Reissner's 
paper,  certain  conventional  simplifications  usually  employed  in  the 
static  loading  case  had  been  introduced  because  of  the  thinness  and 
shallowness  of  the  structure;  the  frequency  equation  was  expressed  in 
the  form  of  a  determinant  involved  with  Bessel  functions;  numerical 
results  were  obtained  by  a  Galerkin  approximation  method.   A  great 
contribution  to  the  problem  of  transverse  vibration  of  thin  shallow 
elastic  shells  was  also  due  to  Reissner  in  1955.   In  his  paper  (46) , 
by  an  order  of  magnitude  analysis,  Reissner  justified  an  important 
simplification  for  the  problem;  i  e  ,  the  tangential  (longitudinal) 
inertia  terms  may  be  omitted  with  negligible  errors.   This  simplifi- 
cation has  made  possible  the  solutions  of  other  shell  oscillation 
problems  (47,  48  and  49) . 

The  work  of  Eric  Reissner  and  others,  as  mentioned  above, 
dealt  mainly  with  linear,  free  vibrations  of  a  thin  shallow  elastic 
spherical  shell   The  first  investigation  of  the  problem  of  forced 
vibration  and  the  problem  of  stability  of  such  a  structure  were 
probably  due  to  Grigoliuk  (11)  ,   In  the  work  by  Grigoliuk,  a  non- 
linear oscillation  system  was  considered  for  the  first  time. 

A  great  amount  of  work  in  the  nonlinear  vibration  of  shell 
structures  and  their  stability  under  periodic  forces  has  been  done 
by  V.  V.  Bolotin.   In  his  paper  of  1958  (4),  he  first  discussed 
the  problem  of  forced  oscillation  and  stability  of  a  complete 
spherical  shell  under  a  periodic  external  loading.   He  considered 
the  oscillations  of  both  nonlinear  and  linear  cases,  that  is,  the 
system  vibrating  with  both  finite  and  infinitesimal  amplitudes  and 


obtained  the  resonance  curves  showing  that  "hard  excitation"  occurs 
until  the  critical  frequency  is  reached.   This  is  in  contrast  with 
the  oscillation  of  a  flat  plate  where  "soft  excitation"  starts  at 
the  lower  boundary  of  the  instability  region.   A  rather  complete 
collection  of  the  problems  of  parametric  stability  of  elastic 
systems  was  also  due  to  Bolotin   His  two  books,  (50)  and  (51) , 
are  no  doubt  valuable  contributions  to  this  field. 

Along  with  the  development  of  aerospace  science,  another  type 
of  shell  buckling  problem  has  become  ever  so  important,  i.e=,  buck- 
ling under  impulsive  and  blast  loadingSo   As  an  extension  of  Hoff's 
work  (14)  on  the  stability  problem  of  a  column,  in  1958,  A.  S.  Vol'mir 
solved  a  problem  of  dynamic  buckling  of  a  hinged  cylindrical  panel 
under  axial  pressure.   In  his  1958  paper,  Vol'mir  considered  the  end 
shortening  of  the  structure  as  a  linear  function  of  time.   Using  a 
first  approximation  of  deflection  form  for  both  total  and  initial 
deflections,  plus  a  Galerkin  method,  he  was  able  to  solve  the  problem 
by  considering  a  1-degree-of-f reedom  nonlinear  oscillatory  system. 
Results  were  obtained  by  numerical  integrations  and  the  critical  condi- 
tion of  the  system  was  determined  from  the  response  curves.   As  an 
example,  he  showed  that  the  dynamic  critical  load  was  about  1.7  times 
higher  than  the  upper  static  critical  load.   The  problem  of  buckling 
of  a  cylinder  under  external  uniformly  distributed  load  increased 
linearly  with  time  was  considered  by  V.  L.  Agamirov  and  A.  S.  Vol'mir 
(V)    in  1959,  using  a  similar  technique  as  in  the  above  mentioned 
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paper  by  Vol'mir   This  same  problem  was  considered  again  by 
Kadashevich  and  Pertseve  (18)  in  1960.   In  contrast  to  Agarairov  and 
Vol'mir's  work,  where  transverse  inertia  was  considered  alone,  these 
two  authors  have  also  considered  the  inertia  of  axisymmetric  com- 
pressing of  the  shell.   By  removing  the  restriction  on  radial  dis- 
placement presented  in  (V)  ,    they  considered  a  nonlinear  dynamic  system 
of  three  degrees  of  freedom   Three  types  of  dynamic  load  were 
considered  in  this  paper:  a  suddenly  applied  load  with  constant 
pressure,  an  impulsive  load  with  finite  time  durations,  and  a  uni- 
form load  increasing  linearly  with  time.   For  a  very  rapidly  applied 
dynamic  loading,  they  found  the  contribution  of  the  inertia  of  axi- 
symmetric compression  is  essential. 

In  1962,  two  papers  treating  the  problem  of  dynamic  buckling 
of  shallow  spherical  shells  under  uniformly  applied  impulsive  load- 
ings appeared  in  the  open  literature  in  the  western  world.   The  first 
paper  was  given  by  J.  S  Humphreys  and  S.  R  Bodner  (15) ,  where  the 
critical  condition  of  the  system  was  determined  by  an  energy  method, 
and  nonlinear  strain-displacement  relations  were  employed   To  the 
author's  knowledge,  this  paper  is  unique  in  that  it  presents  a  solu- 
tion of  dynamic  buckling  problems  in  this  fashion.   Although  the 
behavior  of  the  dynamic  response  could  not  be  obtained  by  this  method, 
the  general  relation  between  the  geometric  shape  and  critical  impulses 
as  well  as  critical  deflections  were  obtained  through  a  much  easier 
and  clearer  discussion   The  second  paper  was  published  by 
Be  Budiansky  and  R.  S.  Roth  (,8)    in  December  1962.   Snapping  of  a 


shallow  spherical  shell  under  an  impulsive  loading  with  various 
time  durations  was  used  as  an  example  to  establish  a  buckling  cri- 
terion proposed  by  these  two  authors.   Using  a  higher  order  of  sym- 
metrical modes,  Budiansky  and  Roth  integrated  a  five-degree  of 
freedom  dynamic  system  numerically,  and  the  critical  condition  of 
the  system  was  determined  by  a  proposed  buckling  criterion  which 
has  a  physically  significant  basis.   It  should  be  noted  that  before 
Budiansky  and  Roth,  the  condition  of  the  threshold  of  the  shell 
buckling  was  determined  in  a  rather  arbitrary  fashion;  i.e.,  having 
neither  physical  nor  mathematical  reasonings.   The  future  trend  of 
the  investigation  of  the  buckling  behavior  of  a  dynamically  loaded 
shell  seems  to  be  such  that  the  governing  differential  equations 
are  Integrated  by  various  numerical  means.   It  is  for  this  reason 
that  the  contribution  of  Budiansky  and  Roth  is  of  particular  interest. 

2,   The  Scope  of  the  Present  Research 

It  is  well  known  that  the  investigation  of  dynamic  buckling 
of  a  shell  under  impulsive  or  blast  pressures  is  mainly  on  determining 
the  magnitude  of  buckling  deflection  in  the  process  of  loading. 
Usually,  on  one  hand,  we  have  a  dynamic  system  which  is  essentially 
"unstable";  i.e.,  the  deformation  is  unbounded  as  time  increases,  and 
we  are  at  the  position  to  determine  the  threshold  of  the  buckling 
from  the  response  curves;  on  the  other  hand,  we  have  a  dynamically 
stable  system,  and  we  have  to  determine  whether  the  buckling  of  the 


structure  occurs  or  not,  because  of  the  significantly  larger  displace- 
ment In  the  transient  region.   The  conventional  method  where  the 
critical  condition  was  determined  by  the  first  amplitude  of  the 
response  curve,  or  by  the  state  where  the  dynamic  response  has  severe 
changes,  may  be  used  to  solve  the  first  type  of  problem  although  it 
will  not  be  precise.   It  may  introduce  serious  errors  if  the  problem 
is  of  the  nature  of  the  second  type   As  we  have  mentioned,  a  buckling 
criterion  has  been  proposed  by  Budiansky  and  Roth  for  a  certain  type 
of  dynamic  load.   These  two  authors  have  based  their  work  upon  a 
certain  physical  picture  of  the  deformation  of  the  shell  during  the 
loading  process  and  established  a  certain  "measure"  which  defines  the 
critical  state  of  the  structure  in  a  characteristic  load* response 
diagramo   The  philosophy  of  this  criterion  is  fresh  and  remarkable; 
however,  the  difficulties  in  extending  to  the  other  cases  is  obvious 
In  a  rigorous  manner,  the  correct  "measure"  of  buckling  can  only  be 
obtained  from  the  sample  of  experiments.   Errors  introduced  from 
plausible  assumptions  may  become  significant  when  the  load- response 
curve  does  not  present  a  change  in  the  form  of  a  "jump."  Therefore, 
the  following  questions  arose:   Is  it  possible  to  propose  a  buckling 
criterion  with  a  more  general  sense?  Is  there  any  relation  at  all 
between  the  dynamic  load-deflection  curve  and  its  counterpart  in  the 
static  case?  Let  us  seek  positive  answers  to  both  of  these  two 
questions   It  is  on  this  basis  that  the  present  research  will  be 
devoted  to  the  following  two  purposes: 


A„   A  dynamic  buckling  criterion  in  a  general  sense;  i.e., 
independ>_nt  of  the  type  of  loading  and  geometry  of  the  shell,  will 
be  proposed  from  a  comparative  basis.   Naturally,  the  critical  con- 
ditions of  the  statically  loaded  shell  and  a  free  oscillating  shell 
vill  be  good  measure  of  the  critical  condition  of  the  shell  under 
dynamic  loadings. 

B.   An  attempt  will  be  made  to  unify  the  two  methods,  i.e., 
the  energy  method  and  the  dynamic  response  method,  in  the  study  of 
the  dynamic  instability  of  shell  structures.   Therefore,  a  qualita- 
tive discussion  of  the  motion  of  the  dynamic  system  as  referred  to  the 
change  of  total  energy  level  is  desiredo 

Because  of  the  second  purpose  mentioned  above,  we  shall 
restrict  ourselves  to  the  problem  of  considering  a  single  deformation 
mode,  i.e.,  a  first-degree-of-f reedom  system.   It  is  well  known  that 
in  the  static  case,  the  above  restriction  will  make  the  result  of  the 
theory  applicable  only  for  sufficiently  shallow  shells,  e.g.,   X<  6 
where 

A^=  ^I2CI-V^)   "'"Vr.K^  4V3(-v>)   Vh    ^'^ 

which  is  a  standard  geometrical  parameter  used  in  the  shallow  shell 
theory  (16) .   It  should  be  mentioned  that,  in  this  research,  we 
shall  use  a  different  geometrical  parameter,  which  has  the  following 
definition: 

e  =   '•7z„  .  (2) 


■Refer   to   Figure    I, 


By  using  this  parameter,  we  shall  have  our  dynamic  equations  in  the 
simplest  form,  i  e.,  ^  appears  only  in  the  linear  term  in  the 
differential  equation.   However,  as  shown  in  (2),   this  parameter 
is  rather  ambiguous.   Therefore,  for  a  proper  interpretation,  we 
always  consider  Q    related  with  X    by  the  following  equation 

e  «   4.V3(j-i/')    /  x^     ,  (3) 

or,  for  V  —     0-3    , 

e    ~   G-609  /  X^  (^) 

hence,  a  larger  Q    value  implies  a  shallower  shell.   It  will  be 
seen  later  that,  for  non-shallow  shells,  Q    becomes  an  awkward 
measure  for  the  geometrical  shape.   In  the  region  of  the  validity 
of  the  present  theory,  however,  the  9    defined  in  equation  (2)  may 
be  satisfactorily  used  as  a  geometrical  parameter.   The  transitional 
instability  of  a  shallow  spherical  cap  under  three  types  of  impul- 
sive loadings  will  be  investigated.   The  dynamic  equations  of  the 
system  will  be  obtained  by  using  a  Galerkin's  approximation  method, 
and  a  phase  plane  method  will  be  employed  to  discuss  the  stability 
of  the  system.   However,  for  the  case  that  the  load  is  a  linear 


^Numbers  in  parentheses  in  the  text  which  are  not  underlined 
refer  to  the  equation  numbers. 


function  of  time,  numerical  integrations  will  be  used  because  the 
system  is  then  nonautononous . 

An  extension  of  the  same  technique  to  the  solution  of 
problems  in  the  case  of  a  complete  spherical  shell  will  also  be 
made. 


CHAPTER  I 

A  CRITERION  FOR  DYNAMIC  BUCKLING 

As  well  as  in  the  static  case,  the  loss  of  stability  of  a 
shell  determined  by  the  defonnatton  of  the  structure,  i.e.,  a 
buckling  state,  can  be  determined  from  the  load-deformation  rela- 
tions.  In  the  static  case,  to  determine  the  state  of  buckling 
usually  does  not  add  any  trouble,  because  there  always  is  a  point 
of  relative  maximum  on  the  load-deformation  curve,  which  is  shown 
by  the  increasing  of  deflection  with  a  decreasing  load,  hence  the 
state  of  instability  is  very  clear   In  the  dynamic  case,  however, 
as  we  know,  the  shell  may  have  buckled  before  the  system  reaches  its 
first  amplitude  in  the  nonlinear  oscillation  due  to  a  severe  change 
in  the  deformation.   Therefore,  a  stable  oscillation  can  cause  the 
critical  condition  for  the  structure  as  well  as  an  unstable  oscil- 
lation.  This  situation  happens  particularly  when  the  load  is 
rapidly  applied, and  with  short  time  durations,  or  an  impulsive 
type.   It  is  for  this  reason  that  a  physically  significant  buckling 
criterion  should  not  be  based  upon  the  stability  nature  of  the  oscil- 
latory system  alone;  it  should  be  safeguarded  by  a  certain  fixed 
value  which  satisfactorily  measures  the  danger  of  the  structure. 


10 


11 


Hovever,  it  is  well  knovn  that  an  unstable  vibration  system  will  have 
its  amplitude  increasing  indefinitely  with  time   Therefore,  the 
transitional  point  for  an  original  stable  system  to  an  unstable  system 
will  always  represent  a  critical  condition  for  the  structure.   This 
is  to  say  that  the  study  of  dynaniic  stability  of  the  oscillatory 
system  is  still  the  most  important  consideration  in  the  investigation 
of  the  dynamic  buckling  of  shell  structures,  although  it  becomes 
impossible  for  some  cases;  for  exanqile,  when  the  system  is  eventually 
unstable,  then  other  techniques  have  to  be  used. 

It  is  the  purpose  of  this  section  to  establish  a  new  buckling 
criterion  based  upon  the  very  nature  of  the  dynamic  stability 
theorems   Certain  measures  of  the  buckling  of  the  shell  of  this 
nature  will  be  provided  after  the  following  discussion.   The  danger 
of  overestimatton  of  the  critical  loadings  will  also  be  safeguarded 
through  the  conqjarison  of  the  characteristic  load-deformation  curves 
for  some  structures  under  other  situations,  whose  stability  nature 
are  well  known. 

The  proof  of  such  a  criterion  is  impossible  at  this  stage, 
yet  its  physical  significance  is  not  difficult  to  observe  and  will 
be  established  through  the  examples  given  in  the  following  chapters. 
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1.   Autonomous  Conservative  System 


The  typical  dynamic  equation  of  such  a  system  is  of  the 
following  form 


-^'  =     Uv,  I)    . 


(1.1) 


vhere  A.  is  a  parameter,  e.g.,  the  load  parameter.   An  equivalent 
form  of  (11)  is  the  two  dimensional  system: 

V  =    5 

i  =  -f  c?,  ^) .  (12) 

It  is  well-known  that  the  discussion  of  the  stability  of  all  the 
possible  motions  described  by  (1.1)  is  essentially  the  same  as 

discussion  the  stability  of  the  motion  in  the  neighborhood  of 

2 
certain  isolated  points,  i  e  ,  the  singular  points,   in  the  phase 

plane  of  the  system.   These  singular  points  are  found  by  the  condi- 
tion that  TJ    and  I  vanish  simultaneously,  i.e.,  from  equation  (1  2), 


^Refer  to  (9),  (^) ,  (^9)  ,  and  particularly  (3),  in  which  a 
beautiful  discussion  of  the  "conservative  system"  has  been  given. 

2 
The  names  critical  point  and  equilibrium  points  are  also 

used. 
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%    =      0 

■f(V,   A)  =    0   . 

(1.3) 
The  first  condition  in  (lo3)  merely  saya  that  the  singular  points 
are  located  on  the  yj   -axis  (where  |=  0  )  •   It  is  the  second  condi- 
tion in  (1.3)  that  determines  the  singular  points  in  the  phase 
plane.   For  a  system  as  (1.1),  we  can  have  only  two  types  of  singu- 
larity, namely,  the  center  and  the  saddle  point.   The  trajectories 
around  a  center  and  around  a  saddle  point  have  a  characteristic 
difference,  and  this  is  shown  in  Figure  2. 

In  Figure  2,  Tj ^     and  ^j^   are  centers;  motion  around  these 
two  points  is  described  by  simple  closed  trajectories,  which  is 
stable  in  character.   The  trajectory  passing  through  the  saddle 
point  is  called  a  separatrix,  which,  less  rigorously  speaking,  is 
the  partition  between  two  motions  with  different  characteristics. 
It  is  also  seen  from  the  same  figure  that  a  trajectory  lies  outside 
of  the  separatrix  and  has  a  higher  energy  level  than  the  one  located 
inside  of  it. 

By  virtue  of. the  above  discussion,  we  may  say  that  the  study 
of  the  stability  of  the  dynamic  system  is  essentially  equivalent  to 
finding  the  character  of  the  trajectory  of  the  system,  and  the  loss 
of  stability  of  the  system  is  equivalent  to  the  condition  that  the 
system  moves  on  the  separatrix  in  the  phase  plane.   A  further  exam- 
ination of  the  phase  plane  sketch  will  make  it  clear  that  the  sense  of 
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"loss  of  stability"  mentioned  above  has  the  same  nature  as  the  usual 
dynamic  buckling  criterions,  i.e  ,  the  characteristic  deformation 
undergoes  a  severe  change  (increased) .   In  the  phase  plane,  all  points 
where  the  trajectories  intersect  the  ^7  -axis  reflect  the  amplitudes 
of  the  motions,  because  j^  =  0  ,  i.e.,  the  '>](z)   curves  of  the 
motions  have  a  horizontal  tangent  at  that  point.   The  magnitude  of 
the  an^litudes  are  measured  relatively  by  the  length  from  the 
origin  0  of  the  phase  plane.   For  any  motion  moving  on  a  trajectory 
inside  of  the  separatrix,  the  amplitude  increases  gradually  as  the 
total  energy  level  increases,  i  e.,  due  to  the  increasing  of  the 
external  pressureSo   This  is  shown  as  from  OA    to  o  g  .   Once  the 
external  pressure  reaches  the  critical  value  which  causes  the  motion 
on  the  separatrix,  the  amplitude  undergoes  a  characteristic  change. 
It  first  reaches  0^7,  ^^^   then  creeps  to  the  magnitude  equivalent 
to  O C  •   Any  motion  outside  the  separatrix  had  its  amplitude 
larger  than  O   C  ,  e.g.,  o I)         The  severe  change  of  the  amplitude 
during  the  loss  of  stability  becomes  apparent  by  comparing  the 
length  of  o  8  with  OC   < 

Let  us  summarize  the  above  discussion  and  make  a  useful 
conclusion.   We  have  reached  the  point  that  the  determination  of 
the  dynamic  instability  of  the  system  (1.1)  is  equivalent  to 
finding  the  motion  on  the  separatrix  in  the  phase  plane  of  the 
system.   It  will  become  clear  in  the  later  examples  that  the 
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equation  of  the  separatrix  is  determined  solely  by  the  unstable 
singular  point  of  the  phase  plane,  hence,  by  one  of  the  roots  of 
equation  (1.3)   Comparing  (1.3)  with  (1.1),  we  immediately  found, 
by  its  very  nature,  that  equation  (1.3)  is  simply  the  state  of  static 
equilibrium,  i.e.,  the  counterpart  of  equation  (1.1)  in  the  static 
state.   Moreover,  we  have  V,     ,  V^    ^^^d  ^j  ,  the  possible  states 
of  static  equilibrium,  in  the  phase  plane  sketch   By  possible  states 
of  static  equilibrium,  we  mean  the  deformation  (or  deflection) 
determined  by  the  position  of  these  points  would  be  a  state  of 
static  equilibrium  if  the  external  distrubance  is  a  static  one. 
Thus  far,  we  are  able  to  state  that  the  loss  of  dynamic  stability 
is  characterized  by  the  load- deformation  relation  reaching  a 
possible  state  of  static  equilibrium.   In  most  problems  of  dynamic 
buckling  of  shells,  the  singular  points  are  interior  to  a  closed 
path.   There  is  a  theorem  due  to  Poincare; 

In  a  conservative  system,  the  singular  points  interior  to 
a  closed  path  are  saddle  points  and  centers.   Their  total  number  is 
odd  and  the  number  of  centers  exceeds  the  number  of  saddle  points 
by  one. 

By  virtue  of  the  above  theorem,  since  in  most  of  the  cases 
of  the  shell  buckling,  the  first  equilibrium  position  always 
corresponds  to  the  trivial  solution  of  the  undeformed  state, 
we  may  state  a  criterion  for  the  instability  of  the  system  has  a 
nature  as  equation  (1.1),  which  is  as  follows: 
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Criterion.   The  threshold  of  the  dynamic  instability  (or 
buckling)  is  defined  by  a  point  on  the  characteristic  load- 
deformation  curve,  where  the  deformation  of  the  dynamic  system 
reached  the  first  unstable  state  of  static  equilibrium. 

It  is  noted  that,  for  a  single  degree  of  freedom  system 
(1.1),  this  criterion  of  instability  should  give  the  same 
result  as  would  be  obtained  directly  from  the  dynamic  stability 
theorems,  i.e.,  the  phase-plane  methods   However,  there  is  no 
restriction  in  the  application  of  the  above  criterion  to  the 
systems  of  higher  degrees  of  freedom,  while  the  topological 
method,  in  general,  does  not  apply  in  such  cases. 

2.   Nonautonomous  System 

In  general,  the  topologic  method  cannot  be  used  to  solve 
the  problem  of  a  nonautonomous  system,  i  e.,  when  the  time  variable 
X  expressly  appears  in  the  dynamic  equation,  because  the  trajec- 
tory of  a  motion  is  in  a  space  rather  than  in  a  plane.   For  a 
certain  class  of  equations,  Minorsky  (24)  developed  a  method  which 
he  called  the  "stroboscopic  method."  By  finding  an  identical  trans- 
formation, the  original  nonautonomous  system  can  be  transformed  into 
a  stroboscopic  system  which  is  autonomous.   Therefore,  the  stability 
problem  of  a  periodic  motion  of  the  original  system  is  equivalent 
to  the  problem  of  investigating  the  stability  of  singular  points  in 
its  stroboscopic  system.   Unfortunately,  this  clever  method  cannot 
be  applied  to  the  type  of  problem  which  has  nonperiodic  motions  and 
with  large  nonlinearity,  mainly  due  to  the  difficulties  of  finding 
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the  stroboscopic  transformation.   Furthermore,  for  certain  problems, 
in  which  we  are  interested,  the  motion  is  known  to  be  unstable  as 
time  increases  indefinitely.   As  far  as  buckling  is  concerned,  we 
are  merely  interested  in  knowing  where  the  deformation  begins  to 
increase  violently  or  attains  dangerous  magnitude.   We  have  seen 
in  the  last  case,  i.e.,  the  autonomous  system,  that  the  beginning 
of  the  violent  increment  of  deformation  is  defined  by  the  initia- 
tion of  the  instability  of  the  dynamic  system,  and  as  a  matter  of 
fact,  they  are  identical.   However,  it  is  impossible  to  'xtend  the 
same  logic  to  the  nonautonomous  system,  for  some  of  them  eventually 
reach  a  state  of  "unstable  motion,"  e.g.,  a  system  under  a  forcing 
function  with  a  magnitude  increasing  linearly  with  time.   No  matter 
what  conditions  we  have,  however,  the  projection  of  the  trajectory 
of  motion  onto  a  ^  -  ^7  plane  still  offers  us  some  information 
regarding  the  "violent  increment  of  deformations,"  as  we  shall  see 
in  the  following„ 

Let  us  take  the  following  system 

TJ    =         -f  (  7]   ,   Z  )      ,  (1.4) 

where  X   is  the  time  variable.   The  similarity  between  (1.4)  and 
(11)  is  easily  obtained  by  taking  T   equal  to  some  definite  value 
of  time,  say  \    ;  i.e.,  at  the  certain  time  X   ,  the  motion  of  (1^4) 
is  on  a  trajectory  in  the  phase  plane  characterized  by  (Id). 
Therefore,  the  motion  of  (1.4)  can  be  treated  instantaneously  as 
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a  motion  of  an  autonomous  system  of  the  form  of  (1.1).   However, 
it  passes  only  one  point  on  the  trajectory  of  each  phase  plane. 
Let  ue  specify,  furthermore,  that  the  form  of  f (  ^ ,  T)  in  (1.4)  is 
increasing  in  magnitude  together  with  t,  i.e.,  the  energy  level 
becomes  higher  and  higher  as  T  increases.   The  space  trajectory  of 
the  motion  of  (1.4),  in  this  case,  can  be  visualized  as  in 
Figure  3.   We  can  project  the  trajectory  onto  a  plane  similar  to 
the  phase  plane  and  it  will  be  in  a  form  as  shown  in  the  above 
mentioned  figure.   It  should  be  remembered  that  we  have  specified 
the  forcing  function  to  be  a  monotonous  increasing  function  of  f. 
Notice  the  form  of  the  trajectory  in  the  "^-"7  plane.   It  is  very 
similar  to  the  form  of  an  autonomous  system  with  a  negative  damping 
terra;  with  the  only  difference  that  the  "unstable  focus"  changes 
with  time.   All  the  unstable  singular  points  ^7^  at  time  tr  =  T^ 
are  determined  by 

■f  ( yj  ,-^i)  =      0      .  (1-5) 

this  is  shown  in  the  figure  as  "^  ,  ^  t    '    '    '   V-  ■      By  virtue  of 
(1.5),  the  following  equation, 


■f  (y)  ,z)    =      0     , 


(1.6) 


is  simply  the  locus  of  all  singular  points  (both  stable  and  unstable 
singular  points)  in  77- J^  plane.   The  intersection  point  of  the  curve 


We  prefer  not  to  use  the  term  "phase  plane"  in  this  case. 
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defined  by  equation  (1.6)  and  the  response  curve,  i.e.,  the  solution 
of  equation  (1.4)  is  simply  the  inflection  point  of  the  response 
curveo   At  the  first  inflection  point  of  the  curve,  the  change  of 
the  slope  is  zero  and  the  slope  of  the  curve  is  a  raaxiinum;  there- 
fore, it  is  the  upper  measure  of  the  "violent  increment  of  the 
deformation."  There  is  another  significant  singular  point  from 
(1.5),  i.e^,  when  t^  =  r  =  0,  the  initial  time.   The  first  unstable 
singular  point  in  the  plane  when  t  =  t^  certainly  is  the  lower 
measure  of  the  critical  deformation^   In  the  case  that  t  appears 
only  in  the  forcing  function,  the  lower  measure  in  the  above  defined 
sense  is  simply  the  critical  amplitude  of  the  free  vibration  of  the 
system  (1.4).   Therefore,  the  middle  point  of  these  two  bounds  is  a 
reasonable  measure  of  the  dynamic  buckling. 

As  a  conclusion,  we  summarize  the  criterion  proposed  in  the 
last  section  in  the  graphs  shown  in  Figure  4.   In  Figure  4  P  repre- 
sents characteristic  load  and  S    is  designated  to  be  the  characteristic 
deformation.   D  represents  the  typical  dynamic  curves,  while  S  repre- 
sents typical  static  curves;  A|_  is  the  critical  amplitude  of  the 
free  oscillation  of  the  dynamic  system.   The  critical  condition  of 
each  system  is  determined  by  Aqy   »  or  bounded  by   A^  and   A^t. 
according  to  the  above  discussions. 


CHAPTER  II 


BUCKLING  OF  A  CLAMPED  SHALLOW  SPHERICAL  SHELL 
UNDER  A  PURE  IMPULSE 


1 .  A  Qualitative  Discussion  of  the  Loss  of 
Stability  of  the  Stiucturc 


It  has  been  shovm  in  Appendix  II  that  the  dynamic  problem 
of  a  clamped  shallow  spherical  shell  under  the  action  of  a 
uniformly  distributed  load  q(t)  can  be  reduced  to  a  single  non- 
linear second  order  ordinary  differential  equation  if  a  first 
order  approximation  of  deflections  in  the  following  form  is 
employed;  thus 


or 


(t)  =  <rt) 


(2„0) 


This   second  order  ordinary  nonlinear  differential   equation,      when 

2 
l)  =   0.3,    has   the   following   form: 


12.1  h^    _^      21^  \       _     10. 6Z        2 


RoV 


53    J_ 
/6    pK 


(2ol) 


The  derivation  of  this  equation  is  attached  in  Appendix  II. 
-Refer  to  Figure  1  for  the  definitions  of  each  quantity 


in  this  equation. 
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If  the  shell  is  sufficiently  shallow  so  that  the  middle  surface  of 
the  shell  can  be  approximated  by  a  paraboloid,  thus 

2  ^  Zo  (  r  /  ro  )  2 

and  the  curvature  of  the  middle  surface  has  to  satisfy  the 
condition: 

—  =   1 
These  last  two  equations  give  us  the  useful  relation: 


Zo  -   ^^7 


2R, 


(2.2) 


By  virtue  of  equation   (2.2),    a  diraensionless   form  of 
(2.1)    can  be  written  as: 


+  2(075  y]^ 

where  the  nondimensional  quantities: 


fn, 


^=\rlp  ^/«o  . 


(2.3) 


h 


Zo   . 


R 


in)=   c-Tz  U*^ 


Eh"-       ^'~'  (2. A) 

have  been  used,  and  P   is  the  mass  density  of  the  shell  material; 
the  dot  (•)  now  represents  d/dt. 

In  that  case  the  external  pressure  is  an  ideal  rectangular 
impulse  as  shown  in  Figure  6.   For  sufficiently  small  ^ T ,  as 
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discussed  in  reference  (28) ,  this  problem  can  be  formulated  in  the 
following  fashion:  the  oscillation  of  the  system  is  essentially 
a  free  oscillation  started  at  the  time  t  =  AT  (i.e.,  A  t -»  O  ),  with 
certain  initial  conditions  derived  from  the  external  disturbance. 
It  is  conventional  to  take  the  displacement  at  the  new  initial 
time,  i.e.,  t  =  ^  T,  still  being  the  same  as  the  one  at  t  =  0. 
For  an  initially  undisturbed  system,  the  new  initial  deflection 
will  remain  zero  because  of  the  smallness  of  ^T.      The  velocity  at 
t  =  AT,  however,  will  not  be  the  same  as  that  at  t  =  0,  because 
of  the  sudden  and  instantaneous  external  disturbance.   If  the 
system  is  originally  undisturbed,  i.e.,  velocity  is  zero  at 
t  =  0,  then  the  velocity  at  the  new  initial  time  can  be  easily 
found  by  use  of  "the  principle  of  linear  impulse  and  linear  momen- 
tum," i.e., 


■^o  i  J   • 


It  is  noted  that  equation  (2.3)  is  equivalent  to  a  nonlinear 
spring-mass  system  with  a  unit  mass  (m  =  1)  and  a  constant  force 
in  its  nondimensional  form,  —  $9      •      Therefore,  the  above  equa- 
tion reads: 


35 

i6 


e  J  At    :=    r^C^^)  ~  0  1 


where  AZ     is  the  dimensionless  form  of  the  time  duration  AT,    by 
the  definition  in  (2. A). 
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From  the  above  discussion,  we  have  our  new  initial  condi- 
tions, thus, 
at  t  =  nt:  , 

"li^Z)   -       0       ; 

where  I*  =  ({'(^^)    is  the  dimensionless  impulse. 

Because  of  the  smallness  of  the  time  duration  of  the  impulse 
(  4X  =5;  O  )  ,  it  is  conventional  to  solve  this  problem  by  taking 
(2.5)  to  be  the  initial  conditions  at  tr  s  0,  thus 

7}  (0)  ^  0  ,  1^(0)=^  ei*  , 

(2.6) 
and  a  free  oscillatory  system  with  a  dynamic  equation  obtained 
from  (2.3)  by  dropping  the  terra  on  the  right-hand  side,  which  is 
involved  with  the  external  load  q.   Thus 

7]*  +    [(276  +5  02^9")^^   -  ^-31  rj^  ■^^   2lo7^rj^j   -  q.^^   ^^ 

The  oscillatory  motion  of  an  autonomous  system  characterized  by 
equation  (2.7)  and  initial  conditions  (2.5)  can  be  discussed  quali- 
tatively by  a  phase  plane  method. 

We  shall  follow  reference  (23)  and  write  (2.7)  in  the 
following  form: 

V  =      I 

i     =  -   C2.76  +3  025  9^)7;  -  5-3}  '7'+  :?'o?  '7M  .   (2.8) 
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Here  it  is  apparent  that  ^  is  the  dimensionless  velocity  of 
the  motion.   Equation  (2o8)  can  be  integrated  once  and  becomes 


-  C 


(2.9) 


where  C  is  an  integration  constant.   Equation  (2.9)  is  the  equa- 
tion of  the  trajectories  In  the  phase  plane.   It  describes  all 
possible  motions  of  a  system  characterized  by  (2.8)  „   The  stability 
of  the  motion  will  be  analyzed  in  the  following. 

The  singular  points  (or  critical  points)  of  (2.8)  are  found 
from  the  equation  (2.8)  by  putting  rj    and  ^  equal  to  zero  simul- 
taneously, i.e. 

^    =      0  ,  and 

77  [2/08  y  -  5' 3 1  V  +  (2-76  +  3  02^  e^)  ]  =    0, 


(2.10) 


It  is  clear  that  the  critical  points  are  located  on  the  '7  -axis, 
One  of  them  is  the  origin  of  the  phase  plane  and  the  other  two 
points  are  defined  by  the  equation: 


v  = 


4-2/6 


^.31  ±  ^o5  Vo'95o3  -  e^ 


(2.11) 


The  stability  nature  of  these  singular  points  can  be  determined 

1 
from  the  following  characteristic  equation,  (9),   i.e., 


Refer  to  p.  317,  reference  (9) 
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■X  ) 


Det  (A)   = 


-[(2.76  +  3  02^  e^)-  10  62^^ 

+  6-?24  X'  ]  _^ 

=        O  . 


(2.12) 


It  is  given  by  the  stability  theorem  that  corresponding  to  a 


pair  of  real  eigen  values  of  (2.12),  an  unstable  singular  point 
will  satisfy  the  following  condition: 

[  6-324  l'    -   10.62  ^,    +(276+3  02^  e')]    <    0.  (2.13) 

It  is  also  apparent  from  equation  (2.12)  that  there  are  only  two 
types  of  singular  points  for  this  system,  namely,  the  stable 
critical  points  of  "center"  type  and  the  unstable  points  of  the 
type  of  "saddle  points." 

Let  us  return  to  equation  (2.11).   It  is  clear  that  there 
exist  real  positive  nontrival  values  of  ^  (singular  points) 
provided, 

G^    6:     0-19305  ,        or 

Physical  significance  of  this  condition  is  that  if  a  shell  is 
sufficiently  shallow  such  that  its  geometrical  parameter  6  is 
larger  than  a  certain  limit,  i.e. 


This  is  the  characteristic  equation  of  the  linear  approx- 
imation.  The  discussion  of  using  this  approximation  is  referred  to  in 
Appendix  lo 
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e   >   0-44,       1  (2.15) 

There  will  be  no  "snap  buckling"  under  the  action  of  an  impulseo 
For  this  reason,  we  shall  be  interested  only  in  those  shells  with 
geometrical  parameter   0  i  O-U-i^,      By  virtue  of  (2.13),  we  know  the 
root: 

(2.16) 
of  (2.10)  is  the  unstable  saddle  point.   The  trajectory  passing 
through  this  point  is  called  a  separatrix.   The  motion  on  the 
separatrix  is  essentially  unstable,  and  the  motion  described  by  a 
trajectory  inside  of  the  separatrix,  in  general,  has  different 
character  than  the  one  described  by  a  trajectory  outside  of  the 
separatrix.   Therefore,  the  problem  of  determining  the  critical 
condition  of  the  system  reduces  to  one  of  find  the  motion  whose 
trajectory  is  the  separatrix. 

Because  the  separatrix  passes  through  the  singular  point 
defined  in  (2.16),  then  by  using  equation  (2.9)  we  can  find  its 
equation.   For  the  purpose  of  emphasis,  let  us  replace  the 
symbol   '^         by   ^   ,  i.e.,  from  now  on,   ?7 '    in  equation  (2.16) 
is  read 

ri'Jl  — 


'or 


This  condition  is  equivalent  to  A  ^387,  where  A  is 


defined  in  equation  (A) 
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From  (2.9),  for  the  separatrix,  we  have 

I  0^38  ^;;  -  3M  \l     +  (2.76  +  3-02^  e')  rj^l 

-^    ]      =       0.  ^2.17) 

It  must   be   remembered   that    from   (2,10),     '^^^also   satisfies   the 
condition: 

(276  +   3  02&  9^)      =  ^.31  \^  -    2loS  -nt     . 

Using  this  equality  and  (2,17),  we  can  determine  the  constant 
C  =  Cs,  which  will  yield  equation  (2o9)  as  the  equation  of  the 
separatrix.   Thus, 

Cs   =    f.77  r]l    ~     /.o^3g  ^,t    1  (2.18) 

where,  from  (2.16), 

\    =        1-2^94S  -  f  (97^1  V  0-J9303  -  0^   .       (2.19) 
Let  us  return  to  (2.5),  the  initial  conditions  which 
define  the  motion  of  the  system  under  the  action  of  an  impulse. 
From  the  first  condition  in  (2.5),  i.e., 

'^(o)      =       0   . 


It  should  be  noted  that   ^jr  is  subjected  to  the  condition  that 
Tjio.i.e.,    ij         is  nontrivial. 
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It  is  apparent  from  this  condition  that  the  motion  always  starts  at 
the  point  where  its  trajectory  intersects  the  t,    axis.   In  case 
the  motion  is  on  the  separatrix,  i.e.,  at  the  critical  condition, 
from  equation  (2.9),  (2.18),  and  with  condition  (2.5),  we  have  the 
follouing  result: 

(  "/,.  C  of   =  ^'(0)  =  Cs  , 

or 

=    i'^/3ze)   /  1.-]-]  r)l  -  ro&38-  ^,^.     , 

(2.20) 
where  '^^^   is  defined  in  (2.19).   These  two  equations  will  give  the 
critical  impulse  for  any  shallow  spherical  shell  whose  geometrical 
parameter  9  is  known  and  satisfies  the  condition  in  (2.14) „ 

An  example  is  given  by  taking  6  =  0.26,  ( A  «;  5) .   We  shall 
see,  particularly  in  this  numerical  example,  that  the  result  obtained 
by  using  a  phase  plane  method  will  be  the  same  as  obtained  by  using 
the  buckling  criterion  proposed  in  the  last  chapter.   Furthermore, 
the  same  result  may  be  obtained  if  Budiansky-Roth' s  criterion  and 
techniques  in  (8)  are  employed.   It  should  also  be  mentioned  that 
the  result  reported  in  reference  (8)  is  numerically  more  accurate 
than  that  given  by  equation  (2.20)  because  of  higher  order  approxi- 
mations used  by  those  authors.   However,  the  problem  solved  by 


They  used  a  five- degree-of- freedom  system. 
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this  simple  but  precise  method  will  permit  certain  qualitative 
conclusions  which  could  not  be  obtained,  or  would  cause  much  labor 
in  calculations  if  other  methods  are  employed. 

When         e   -     0-26    ,     (2.5),    and    (2.6)    read: 

yj    =     -  1.101^   V(V-  o.S35o3')(r)-  ).68453) 

rjco)    =    o       ,  TJCo)    -      0-^362^  !'■'. 

(2.21) 

Three   singular  points  on   the   Tj  -axis  are: 
^     —      O       ,  a  center 

y]     =      0-J3S03      ,  a  saddle  point 

r]    =      1  6j4fe3     ,  a  center. 

Let      '7t»-     =   0.83503;    from   (2.20),   we   found 
1^^     =        1.34233    .  (2,22) 

The   equation  of   the   separatrix  is   found  as   follows: 

|^=     -  \-o^i%  rf    +  j^4  V^  -  2-9(>i^^  7^  +  0-^1850. 

This  equation  and  other  phase  plane  trajectories  have  been  plotted 

and  are  shown  in  Figure  7. 

•,:r     ''-- 
In  Figure  7,  when  I  <  1  cr   given  in  (2.22),  trajectories 

of  the  motion  are  closed  curves  around  the  center  point  at  the 
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origin,  e.g.,  curve  1;  when  1  =  l^r    .  the  motion  is  unstable  and 
on  the  separatrix  2;  when   1   >  1  cr   »  the  motion  is  on  a  trajectory 
such  as  3o   It  is  also  easy  to  explain  the  occurrence  of  the 
"buckling"  from  this  figure.   The  points  where  trajectories  intersect 
the  7  -axis  correspond  to  the  situation  that  the  response  curve 
reaches  its  amplitude.   Therefore,  the  ^    coordinate  of  these 
points  (e.g.,  OA)  is  the  measure  of  the  maximum  inward  central 
deflections  of  the  shell.   It  is  clear  that  the  maximum  inward 
central  deflection  increases  with  i  in  a  continuous  fashion  when 
1  <  1  cr    •   As  soon  as   I   is  slightly  larger  than   I  <.r   ,  the 
maximum  central  deflection  undergoes  a  severe  change,  from  some 
value  less  than  0.84  (e.g.,  OA)  to  some  value  greater  than  2 
(e.g.,  OC) o   Because  of  this  severe  change  of  deflection,  snap- 
buckling  of  the  shell  occurs. 

2.  A  Study  of  the  Dynamic  Response 

It  is  noted  that  the  differential  equation  has  the  form 
of  (2.7),  and  with  initial  conditions  (2.6)  can  be  integrated.   The 
solution  of  Such  an  equation,  in  general,  is  involved  with  Jacobian 
Elliptical  Functions,   It  is  still  impossible  to  give  a  nontrivial 
expression  for  the  solution  of  the  equation,  which  is  of  the  same 
form  as  equation  (2o7)o   However,  when  9  is  taken  to  be  a  definite 
value,  the  solution  of  (2.7)  can  always  be  obtained.   For  an  example, 
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In  solving  such  type  equations,  we  shall  take  the  system  defined 
in  equation  (2.21);  i.e.,  6  is  taken  to  be  0„26  in  equation  (2.7) 
and  (2.6).   All  numerical  work  involved  will  be  presented  in  detail. 
We  feel  that  the  result  of  this  section  will  clarify  certain  impor- 
tant points  in  both  the  last  section  and  the  following  section  on 
the  justification  of  the  buckling  criterion. 

Let  6  =  0o26  in  equation  (2.9);  we  have  the  equation  of 
the  trajectories  for  the  system  (2.21)  in  the  following  form: 

^  (2.23) 

where  C  is  an  integration  constant.   By  using  the  second  initial 
condition  in  (2.21),  C  has  the  following  expression: 

C  =  (0-5362^   l")^  .  (2.24) 

Our  problem  is  to  find  the  response  '^(t)  corresponding  to  each 
disturbance  I*.   It  is  still  inq)08sible  to  obtain  the  general 
expression  and  only  particular  cases  will  be  given. 

We  shall  study  the  responses  corresponding  to  two  individual 
disturbances:   1  =  1.2,  and  i  =  1.5.   We  shall  see,  in  contrast 
to  the  small  increment  in  the  disturbance  I  ,  the  corresponding 
responses  will  undergo  characteristic  changes. 

When  I  =1.2,  the  positive  branch  of  (2.23)  reads: 

-    =    > 


(2.25) 
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This  can  be  transformed  into  an  elliptical  integral  of  the  Legendre's 
standard  form.    Let  us,  first,  formulate  the  quadratic  equation: 

2.8'6166  V  -  ^iSi(>i>  V    -f  >  /0293  =   0  .    (2.26) 
Its  coefficients  are  related  to  the  zeros  of  the  algebraic  equa- 
tion under  the  radical  sign  in  (2.25)  in  definite  ways  which  can  be 

found  in  almost  any  textbook  treating  on  elliptical  functions  and 

2 
will  not  be  given  here.   The  two  roots  of  (2o26)  are: 

P    =  027585    ,     I-    =       hi9li»    .  ^^'^''^ 

Now,  we  use  the  following  transformation: 


^  = 


1+  Z 
—         o  27Sgr^  ->-  \-59i]2;  Z 

1+2  »  (2.28) 

and 

and  then  equation  (2.25)  can  be  transformed  into  the  following  form: 


1 
Refer  to  reference  (41) 

2 
eg.,  reference  (43) , 
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I   12135       dZ 
(1  +  2)^       ^^ 


1      /l  O^SSfo  y7721  Z^-  i  ^3^GZ)(  1i^2lZZ^-OI69S 3) 


(i  +  2) 


dz 

046039     dz   = 


(2.29) 

Referring   to   reference    (41)/   the   function   z   can  be  written  in  the 

form: 

Z  =      0  34^95       Sn   fu  Jm)  ,  .^  ^o) 

where 

m  =     o-o]u63  . 

(2.31) 

Therefore,  as  we  substitute  (2.30)  into  equation  (2.28),  the 
solution  of  the  problem  can  be  formally  e3q>re3Sed  as: 


1 
The  expressions  are  on  p,  26,  reference  (41) 
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7|  _        o-27fe»5      +     0  4SI96    Sn(^\^) 

I       +      0  34495   S>i(u|7n) 

(2o32) 

where  U  and  m  have  their  definitions  in  (2.31) o   The  value  t 

o 

in  (2.31)  is  determined  by  the  initial  condition:  when  r  =  0, 
y]     =0.   In  this  fashion,  we  have  Zo     satisfying  the  following 
equation, 

Sn  (-130995  Zo  I  0  01469)   =   0-57235. 


As  an  approximation, 

Z.    =      -  0  469^ 


(2.33) 


(2.34) 


Solution  (2.32)  has  the  following  general  properties: 

Ao   It  is  periodic  because  it  involves  the  double  periodic 

function  5n(uliTi')  .   The  real  period  of  SuCulm)  is 

AK  =  4(1.57658).   The  period,  P  of  •»!  ,  then,  is  equal  to  4.81419 

according  to  equation  (2.31). 

B.   The  ^  values  are  bounded  in  the  interval  -0.31464  ^ 

V    -      0o56344,  because  the  value  of  S-n  (u  |m)  varies  between 

-1  and  +1.   Therefore,  the  maximum  amplitude  of  the  dynamic  response 

t)    corresponding  to  the  disturbance  I*  =  1.2  is  ri    =  0.56344. 
'  '  max 

The  next  exaiqjle;  when  I*  =  1.5,  the  positive  branch  of 
(2.23)  reads: 


dr; 


-^  =  ^  =  y  0-64702  -  2964^  ri^   +  3^4  ^^  -  /  o^Sg"  ^' 
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Different  from  the  last  case,  the  rational  function  in  the  radical 
sign  in  the  above  equation  has  complex  roots.   As  a  counterpart  of 
equation  (2.25),  we  can  express  the  above  equation  in  the  following 
form;  thus 


dr) 


0  7/  I 

—  =         J -]o^3S    (7]-ol)  i7]-^)(7]-  r)(   V-  ^)     .    (2.35) 
where 

d  =   2.  o(>^S^      ,  a    =     -   o.3go67  , 

r     =      o  gjt^6  +   0  28'5j7  t     ,  5    =    ogjt^6  -0  28557  i. 

The   counterpart   of   equations    (2.26),    (2.27),  and    (2„28)    are 
respectively 

o  ono3  IT    +    3  I3ij26  V    -2  632/^      =        Oi 

p-     -lSnSl9     »  *&    =       0.8'3&8-7    '  (2.36) 


and 


1      -^     Z 


(2.37) 
The  transformation  in  (2=37)  brings  equation  (2.35)  into  the  simpli- 
fied form  as  in  the  following: 
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where 


A     =     /<o538'(t/j.9993)o.o8iua/j55.7//+^72      , 

M,    =      3/i,49l-  ^60^9  , 
M2    =      34.  ^?9- 97923  . 


m  =    o.94g'49. 

Substituting  (2.39)  into  (2.37),  we  have  the  solution 

-\»k2l9     CnCulm)     +      I26-75JI*I 


1 


where  the  definition: 

Nc  =    V  Cn   , 
has  been  used,  and  u,  m  are  defined  in  (2.40) 


(2.38) 


Therefore,    we  have 

Z    =      \^\GU2i4.^      Ncinim.),  (2.39) 

where 

U.    =      1. 29092    iZ-to)    , 


(2.40) 


C-n(ulm)     +      /5l-64-2i|9  ♦ 

(2.41) 
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In  a  similar  way,  x^  in  (2.40)  can  be  determined  by 
requiring  rj(i=o)=    0     I    thus 

l^     -     -0.70755  (2.42) 

It  is  obvious  that  the  solution  in  equations  (2.41),  (2.40),  and 
(2.42)  is  characteristically  different  from  the  solution  represented 
by  equations  (2.32),  (3.31),  and  (2.34).   Solution  in  (2.41)  has  the 
following  characteristics: 

Ao   It  is  periodic,  with  period  P  =  9.03„ 

Bo   "7  values  are  in  the  interval  -3.8079  i    V    -   2„06868. 

The  maximiim  amplitude  of  the  dynamic  response, 
"Imax  =  2,06868„ 

Notice  the  characteristic  change  in  the  form  of  the  dynamic 

response  and  the  severe  increment  in  the  amplitude  (from 

'V    =  0„563  to  2*068)  as  I*  value  changed  from  1,2  to  1.5.   We  can 
max 

conclude  that  the  critical  load  I^^  must  be  some  value  in  between 
the  two  values.   One  gets  a  satisfactory  justification  by  referring 
back  to  equation  (2.22),  where  the  critical  impulse  was  found  to 
be  1.34233. 

Response  curves  corresponding  to  I*  =  1.2  and  I*  =  1.5  are 
presented  in  Figure  8o 
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3o  A  Justification  of  the  Buckling  Criterion 

In  the  last  section,  we  have  seen  that  the  dynamic  response 
for  a  system  defined  by  (2.21)  can  be  found  by  integrating  the 
differential  equation  directly,  and  the  solution  in  terms  of  Jacobian 
Elliptic  Functions.   The  dynamic  responses  corresponding  to  other 
external  inqjulses  than  those  given  in  the  last  section  may  also  be 
obtained  in  a  similar  manner,  yet  the  procedure  is  laborious.   If 
merely  the  amplitudes  of  the  response  curves  are  desired,  then  for 
a  one-degree-of- freedom  system,  as  equation  (2.21),  the  difficulty 
of  integrating  the  differential  equations  can  be  removed  by  use  of 
the  information  obtained  from  the  previously  discussed  topological 
method. 

We  have  mentioned  that  the  intersection  points  of  the 
trajectories  and  the  V  -axis  in  the  phase  plane  are  the  points 
where  the  response  curves  reach  their  amplitudes,  because  at  those 
points,  the  velocity  ^  is  equal  to  zero.   This  fact  suggests  that 
we  obtain  the  amplitude-impulse  relation  for  the  system  (2.6)  and 
(2.7),  or  their  special  case  (2.21),  from  the  equation  of  trajectories, 
i.e.,  equation  (2.9).   Let  us  restrict  ourselves  to  dealing  with  the 
special  case  defined  in  equation  (2.21);  i.e.,  9  is  taken  to  be  the 
value  0.26.   Substituting  (2.24)  into  equation  (2.9)  and  using 
6   =  0.26,  we  have  the  following  equation  for  the  trajectories 
of  this  system: 
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(2.43) 
By  virtue  of  the  above  discussion,  the  amplitudes  of  the  response 
curves,  ^ijj's,  are  found  by  setting  the  velocities   E,  '  s  equal  to 
zerOo   By  doing  so,  from  (2  43)  the  following  relation  is  obtained: 

(2.44) 
For  the  solution  to  be  physically  meaningful,  the  positive  branch  of 
the  last  expression  should  be  used;  thus. 

This  relation  is  shown  in  Figure  9   Good  agreement  between  the 
result  presented  in  Figure  9  and  the  result  for  the  two  cases  worked 
out  in  the  last  section  indicates  the  correctness  of  this  technique. 

In  Chapter  II  we  have  proposed  that  the  shell  will  buckle 
when  the  characteristic  load-deformation  curve  reaches  the  first 
unstable  state  of  static  equilibrium   There  are  three  particular 
static  equilibrium  states,  i.e,  the  positions  of  rest  of  the  system 
(2.21)  which  can  be  found  from  the  first  equation  of  (2.21)  by 
setting  the  inertia  force,  i  e.,  'tj     term,  equal  to  zero.   Thus,  the 
three  positions  of  rest  (where  ''7  =  "^    ~    %    -  0)    are  obtained: 

11)  12)  '3>      ^  ^ 

rj   -     0  ,  ^    -    0214-     ,      "7  =  /-^g-  .         (2.46) 
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In  '?niax  "  ^*  Pl^'^s*  these  curves  will  be  straight  lines  parallel 
to  the  I*  axis.  According  to  the  criterion,  the  first  unstable 
equilibrium  position  (  ^  =  0.84)  defines  the  critical  conditiono 
As  shown  in  Figure  9,  the  l*^.^   thus  found  is  identical  with  that 
found  previously  by  the  phase  plane  method. 

It  is  also  of  interest  to  see  the  comparison  between  the 
present  criterion  and  that  proposed  by  Budiansky  and  Roth  in  (Q) . 
According  to  Figure  1  and  equation  (2.0),  we  have: 


Therefore,  certain  definitions  in  reference  (8)  assume  the  following 
expressions: 


^2  , 


Jq   J  o 
By  definition  we  have 

A  =       ^/Zo      =  V    .  (2.47) 


According  to  (2.47),  the  measure  of  buckling  used  in  (,8), 

/^  =1,  corresponds  to  17    =  1  in  the  notation  of  this  paper. 

max  max 
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It  is  apparent  from  Figure  9  that  the  same  critical  impulse  will 
be  obtained  if  the  criterion  proposed  in  reference  (8)    is  employed. 
Furthermore,  it  is  also  indicated  in  the  same  figure  that  the  measure 
used  by  Budiansky  and  Roth  falls  into  the  unstable  branch  of  the 
load- deformation  curve  and  is  close  to  the  point  of  instability. 
Therefore,  the  criterion  proposed  by  these  two  authors  is  proved  to 
give  satisfactory  accuracy  for  this  specific  problem. 

Critical  impulses  for  other  values  of  6  based  upon  equa- 
tions (2.19)  and  (2.20)  have  also  been  calculated.   Results  are 
presented  in  Figure  10.   It  is  shown  in  the  figure,  for 
0,44 >  e  >  0.32  (or  3=87  <  X    c    4  53)  that  equation  (2.20)  agrees  well 
with  the  result  given  in  (8)    and  appears  almost  the  same  as  the 
result  of  reference  (15)  ,  when  9<  0.15  (or  A  >  6.6).   It  is  believed 
that  this  analysis  is  parallel  to  the  presentation  of  (15) ,  yet  with 
a  cut-off  point  at  a  larger  ©  value,  i.e.,  the  present  analysis 
admitted  a  shallower  shell  to  buckle  under  the  applied  impulse.   This 
tendency  seems  to  be  correct  as  compared  with  the  result  of  using  a 
higher  degree  of  approximation  given  in  reference  {8)  . 

4.   A  Note  on  the  Effect  of  Initial  Geometrical  Imperfections 

It  is  rather  interesting  that  we  may  conclude,  on  the  basis 
of  Figure  7,  that  any  axially  symmetric  geometrical  imperfection  will 
give  a  deduction  of  the  critical  impulse  for  the  shell.   It  is  the 
nature  of  the  equation  of  the  separatrix  of  having  a  relative  max- 
imum when  "7  ~  0,  i.e.,  on  the  J,  axis.   Any  initial  imperfection 
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(deflection)  of  the  shell  is  equivalent  to  set  the  motion  starting 

at  C  '  instead  of  at  C^  as  shown  in  Figure  7.   It  is  seen  from 
s  s 

this  figure,  that  Cg'  has  a  smaller  ordinate  than  C  .   Therefore, 
the  critical  impulse  based  upon  the  former  will  have  a  lower  value. 
For  an  example,  let  us  assume  that  the  initial  deflection  of  the 
shell  is  axially  symmetrical  and  has  the  same  form  as  the  deflection 
of  the  shell,  i,eo,  can  be  described  by  equation  (2cO);  further- 
more, it  has  a  central  deflection 

or 

'^.  —  00^     _  (2.48) 

Based  upon  this  value,  the  critical  impulse  will  be  6  per  cent  less 
than  that  directly  given  by  equation  (2.20) „ 


CHAPTER  III 


BUCKLING  OF  A  CLAMPED  SHALLOW  SPHERICAL  SHELL  UNDER  A 
SUDDENLY  APPLIED  CONSTANT  PRESSURE  OF  INFINITE 
TIME  DURATION 


1.  A  Qualitative  Study  of  the  Instability  of  the  Structure 

We  have  seen,  from  the  last  section,  the  solution  of  the 
shell  buckling  under  a  pure  impulse  with  infinitesimal  time  dura- 
tion did  not  answer  the  buckling  problem  of  the  shell  under  the 
impulsive  pressures  with  sufficiently  long  time  durations.   In 
this  section  we  shall  solve  the  problem,  which  is  another 
limiting  case  of  the  above-mentioned  problem;  the  time  duration 
of  the  applied  impulsive  loading  is  infinite. 

Let  us  assume,  without  loss  of  generality,  that  before 
the  time  t  =  0  the  system  is  at  rest,  and  at  time  t  =  0  a 
uniformly  distributed  constant  pressure  q  is  suddenly  applied  on 
the  surface  of  the  shell.   The  history  of  such  a  load  is  shown  in 
Part  B  of  Figure  6.   The  dynamic  equation  of  this  problem  can  be 
obtained  directly  from  (2.3)  by  taking  q  as  a  constant o   Thus, 

^   ^        f^H'   [(^76  +  JO250')  7]     -   1.31  ^^ 

+  2108-  7/M      Z  >  0  .  (3.1) 
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We  have  mentioned  that  the  system  is  initially  at  rest,  i.e., 
y)  (O)    =        O 

are  the  conditions  at  t  =  0. 

Note  that  the  autonomous  system  in  (3.1)  can  be  integrated 
once  and  results  in  the  form: 

i,^  -      A  +  (uiz^  9^)  7j  -  (2  76  +  3  02^  e"- )  V^ 

(3.3) 

where  A  is  the  integration  constant.   Equation  (3.3)  is  obviously 
the  equation  of  the  phase  plane  trajectories. 

The  singular  points  (critical  points)  of  the  system  in 
equation  (3.1)  are  found  from: 

'i   =     ^    =     0 

ioe.,    on   the    i]    axis  where 

tI  ef  -    [  (2.76  +5  02^e^)rj   -   531    f 

+  2  lo^  rj'i    1       —  Q_ 

Stability  conditions  in  the  neighborhood  of  these  singular  points 
are  determined  by  the  following  characteristic  equation,  thus 


(3.4) 


;\^  +  [6  72fe5  7]^-  /0  62  r]  +C2  76+  5  0^*5  0^)]     =0, 


3.5) 
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Corresponding  to  the  real  eigen  values  yielded  by  (3  5),  the  saddle 
points  satisfy  the  following  condition: 


6  3225  ^/  -  JO  62  V,    +  (2  16  -^    3  02^    e')     <     0     ,         ,,  ., 

(J.b) 


where  y]      is  the  root  of  equation  (3.4) 

The  problem  of  finding  the  critical  load  q  which  will  make 
the  motion  of  the  system  described  by  (3  1)  and  (3.2)  lose  its 
stability  is  equivalent  to  the  problem  of  finding  the  q  value  which 
will  put  the  motion  of  the  system  on  a  trajectory  passing  through 
the  saddle  point  in  the  phase  plane,  i.e.,  on  the  separatrix.   By 
virtue  of  its  initial  conditions,  the  trajectory  of  the  motion  of 
the  system  (3.1)  and  (3.2)  passes  through  the  phase  plane  origin 
(  7^  =  E,  =  0)  .   Therefore,  the  equation  of  this  trajectory  is 
obtained  from  (3.3)  by  taking  A  =  0,  thus 

1^  =    [  a.i2^  e%  ~   (276  +  3  025  9^)ri 

If,  furthermore,  this  trajectory  also  passes  through  the  saddle 
point  (  y]     ,  0),  it  is  obvious  that  the  following  condition  has 
to  be  satisfied: 

or  equivalently 

f  ef   -   [(?76+  5025  0')^,-  J.^4   V'   +    10>SU    Vo] 

=        ^    '  (3.8) 
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for   "Vo  ^   0. 

Attention  is  invited  to  the  fact  that  this  problem  reduces 
to  finding  common  roots  between  equations  (3.4)  and  (3.8).   A 
solution  to  a  similar  problem  for  a  complete  spherical  shell  is 
given  in  Chapter  V,  where  the  technique  is  discussed  in  detail. 
Equations  (3.4)  and  (3.8)  will  have  common  roots  provided 

_        ITT  yjt    -  lo^u    rj^ 

%      ~  ■  (3.9) 

"  2 062^  e  ' 

Substituting   (3.9)    into   (3.4),   we  have 

3  162    7^'  -   7  08-  77^  +    (  Z76    +  5  02^    9^)   7      =         0, 

(3.10) 

which  has  three  roots.   The  two  nontrivial  roots  are  solved  from 
the  following  equation: 


/  (  /5.2I79Z  -  3^-26020  Q^) 
Tj         =      1  Il95i^  +  / 

'•^        V      3999297 


(3.11) 


It  should  be  noted  that  it  is  not  true  that  both  solutions  in 

(3.11)  are  the  singular  points  of  the  system,  because  they  are 

not  solved  directly  from  (3.4) „   It  merely  says,  at  the  present 


^Let  us  recall  that  V^   in  (3.8)  is  one  of  the  roots 
of  equation  (3.4);  it  satisfies  the  condition  (3.6). 
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time,  that  both  ^7,  and  ^^  in  (3.11)  are  the  possible  conmon  roots 
between  (3.4)  and  (3.8)  vhen  the  corresponding  q  value  in  (3.9) 
is  taken.   In  other  words,  if  'f]^    is  substituted  into  (3.9)  to 
yield  a  particular  value  of  q  and  this  value  is  used  to  replace 
the  parameter  q  in  both  (3  4)  and  (3  8),  then  these  two  equations 
will  possess  a  common  root,  ^7  =   ^,   •   (Note  that  T^p  may  not  be 
a  root  of  either  of  these  two  equations.)   Since  equation  (3.10) 
gives  the  common  root  of  (3  4)  and  (3.8)  which  specifies  the 
buckling  of  the  shell,  we  may  deduce  one  of  the  important  con- 
clusions, i  c  e. , 

d        ^  0-6307  ^,  (3.12) 

has  to  be  satisfied;  otherwise,  no  buckling  will  occur,  because 
equation  (3.10)  has  no  real  positive  nontrivial  roots,  as  indicated 
in  (3.11)« 

We  shall  see  that  only  one  of  the  two  nontrivial  roots  of 
(3.11)  will  satisfy  condition  (3.6),  i.e.,  be  a  singularity  in  the 
form  of  a  saddle  point.   Discarding  the  trivial  solution  Tj  u   0 
from  (3.10),  we  have  the  following  condition: 

{?.76  +302^  0^)  -  7.o»rj-    3  162  rj^    .  (^"13) 


It  corresponds  to  A  >  3.2123,  where  A  is  defined 
in  equation  (4) . 


48 


which  is  satisfied  by  both  fj     and  V    in  (3oll).   Substituting 
(3.13)  into  (3.6),  after  combining  terms  in  the  same  order  of  ^, 
the  following  condition  is  obtained: 

3  162  rj^    -    3^U    Tq       <     0      . 

As  we  mentioned,  both  ^i  and  ^„  are  real  positive  values, 

because  we  have  put  the  restriction  (3.12)  into  our  problem.   If 

this  is  the  case,  and  it  is  noted  that  "H  is  solved  from 

'o 

equation  (3.10),  condition  (3.6)  is  reduced  to  the  following  form: 

y}  <   j.  (19^4  .  (3.14) 

It  is  apparent  that  only  the  smaller  root  in  (3.11)  will  satisfy 
the  condition  (3.14)  and  be  a  possibly  unstable  saddle  point  for 
the  dynamic  system  (3»1).   More  definitely,  let  us  put  the  critical 
deflection  as  follows: 


V        39  99291 


(3.15) 


Therefore,  without  even  going  back  to  the  phase  plane,  we  can 

write  a  general  solution  for  the  critical  pressure  by  substituting 
(3=15)  into  (3.9),  i.e., 

''^  2.  o425  9  .  (3ol6) 
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where  71   has  its  definition  in  (3.15). 
'cr 

It  may  be  apparent  that  certain  techniques  and  their  logical 
foundations  have  not  been  made  clear  in  the  above  discussions,  and 
they  are  very  difficult,  if  not  impossible,  to  be  discussed  on  a 
general  basis.   Therefore,  we  shall  give  an  example  using  a 
specific  shell  with  6  ■  0o26,  which  has  been  used  in  the  previous 
problema 

Let  us  first  find  the  solution,  i.e.,  the  critical  deflec- 
tion and  load,  and  then  go  back  to  verify  that  these  critical 
quantities  do  put  the  motion  of  the  system  on  its  separatrix  and 
cause  the  shell  to  reach  the  threshold  of  instability. 

From  (3ol5),  when  6  =  0.26,  we  have 

and  from  (3.16),  we  obtain  immediately, 

Our  question  has  been  whether  the  value  of  q  in  (3.18)  does  cause 
the  instability  of  the  system  (3  1)  when  Q   =  0.26,  in  other  words, 
to  make  the  motion  of  the  system  on  its  separatrix,  or  equivalently, 
does  the  value  of  T]    in  (3.17)  define  the  saddle  point  in  the  phase 
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plane  when  q  takes  the  value  in  (3.18)?  To  answer  this  question,  let 
us  write  down  equation  (3.4)  and  (3.8)  by  taking  6  =  0  26  and 
q  =  q^-j.  defined  in  (3.18);  thus,  equation  (3.4)  is  in  the  form: 

2.IO?  Tj'^   -  5.?l  rj^   +  2  9(>uk9  ^  ~   0367^3  =  0  ,    (3.19) 
which  has  three  real  positive  roots: 


7j"=  o /75o5  ,    V^=  0  5575U  ,     "7  =  17863?.  (3. 20) 


These  are  the  singular  points  in  the  phase  plane,  while  equation 

2 
(3 08)  has  the  following  form: 

[  I  o5a  y]^  -  3^k  "^  +2  9(,u.^  t]    -  073507]     =      o,     (3.21) 

which  has  three  roots: 

and  a  double  root  at    rj  =    0  i:^7fe^  .      (3.22) 

Comparing  (3.20)  and  (3.22),  it  is  clear  that  "^  ;  0.55754,  i.e., 
the  critical  value  defined  in  (3.17)  is  the  common  root  between 


We  have  multiplied  the  value  (-1)  through  the  original 
equation  (3.4) . 

2 
Similarly,  a  quantity  (-1)  has  been  multiplied  through 

the  original  equation. 
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(3,19)  and  (3.21).   Furthermore,  tested  by  condition  (3.6),  this 
value  (  ''J  =  0.55754)  defines  a  saddle  point  in  the  phase  plane. 
A  further  remark  on  the  equations  (3.21)  or  (3.8)  will  con^letely 
answer  the  above  question.   Let  us  replace  the  Tj     in  the  first 
equation  of  (3.8)  by  "^  ,  and  comparing  the  final  equation  with 
equation  (3,7),  we  find  (3.8)  is  merely  the  condition  of  the  inter- 
section of  a  phase  plane  trajectory  with  the  ''^-axis  (  ?>  =  0)  . 
Because  of  the  preceding  discussion,  we  conclude  that  the  phase 
plane  trajectory  for  the  motion  of  this  particular  shell  (   9  =   0.26) 
will  pass  through  the  saddle  point  if  the  parameter  q  is  defined  in 
(3.18);  therefore,  the  value  of  q  in  (3.18),  or  in  general  in  (3.16), 
is  the  critical  value  for  the  external  pressure.   The  phase  plane 
trajectory  of  the  motion  of  the  shell  at  the  threshold  of  the 
instability  is  shown  in  Figure  11 o   Based  upon  equation  (3.16)  and 
(3.15),  relations  between  the  critical  load  and  the  geometrical 
parameter  Q    are  found  and  given  in  Figure  12.   Curves  indicating 
the  variation  of  critical  deflections  with  0   values  are  given 
in  Figure  13.   Comparison  between  this  problem  and  its  counterpart 
in  the  static  case  has  been  made  and  is  also  shown  in  the 
figures  mentionedo 
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2.   Another  Justification  of  the  Buckling  Criterion 

In  Section  3,  Chapter  II,  we  have  discussed  a  method  of 
obtaining  the  load-deformation  curve.   For  a  different  problem 
in  this  chapter,  the  same  method  may  still  be  used  to  obtain  the 
relation  between  the  characteristic  load  and  the  characteristic 
deformation.   The  critical  condition  of  this  system,  according 
to  the  discussion  in  Chapter  I,  can  be  determined,  and  the  result 
thus  found  will  be  compared  with  that  obtained  in  the  last 
section  as  a  justification  of  the  proposed  criterion. 

If  the  characteristic  deformation  is  chosen  as  the  ampli- 
tude of  the  central  deflection  following  Section  3,  Chapter  II, 
the  amplitude-pressure  relation  for  system  (3.1)  can  be  obtained 
from  the  equation  of  trajectories,  i.e.,  equation  (.3.7).   By 
equating  ^  to  zero,  we  obtain  the  following  relation  between  the 
nontrivial  amplitude  of  the  response  curves  and  the  corresponding 
djmamic  loadings. 

-       (-?  76  +  J  02^  0^  )  t;  -  3^A  "7^  +  (O^A  ^^ 
d  ~  4.'26  e  •       (3.23) 

For  the  purpose  of  further  discussions,  let  us  take  a 
specific  shell  with  9   =  0.26.   From  (3.23),  we  obtain: 


o      _         2.964U9    r}    -  3^U   ^^     +     I  o^U  ^' 


'^  1012^  .         (3.24) 
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The  static  load-deflection  relation  can  be  easily  found  from 
equation  (3.1)  by  taking  the  inertia  term,  i.e.,  the  '^    term  equal 
to  zerOo   For  6    =   0o26,  we  have 

\     ~  0^Z62^  '  (3.25) 

Equations  (3.24)  and  (3c 25)  are  plotted  in  Figure  14.   The  inter- 
section point  of  these  two  curves  is  at  V^  =   0.55754,  where  the 
dynamic  curve  has  a  relative  maximum.   It  is  obvious  that  the  inter- 
section point  falls  into  the  unstable  branch  of  the  static  curve. 
According  to  Case  1  of  the  buckling  criterion  in  Figure  4,  the 
corresponding  pressure  at  the  point  of  intersection  is  the  critical 
dynamic  load,  thus  q^^j.  =  0.68539.   These  results  are  identical 
with  those  given  by  (3.15)  and  (3.16)  based  upon  the  stability 
theorem. 

The  identical  results  obtained  via  two  different  approaches 
have  established  the  following  facts: 

Ao   In  certain  cases  the  dynamic  criterion  proposed  in 
Chapter  I  is  identical  with  the  stability  theorem. 

B.   Without  given  explanations,  we  have  taken  the  specific 
expression  in  equation  (3.9),  which  led  to  finding  the  solution 
of  the  problem.   This  turns  out  to  be  correct,  since  the  solu- 
tion in  (3.23)  is  entirely  independent  of  (3.9),  yet  the  same 
result  was  yielded. 

A  comparison  between  the  present  theorem  and  reference  (30) 
is  presented  in  both  Figure  12  and  Figure  14,   It  is  interesting 
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that  the  critical  dynamic  load  given  by  the  present  theorem  falls 
in  between  the  two  values  for  the  critical  dynamic  load  obtained 
by  use  of  different  methods  given  in  reference  (30) .   The  critical 
deflections  given  by  (30)  are  seen  to  be  larger  than  those  given 
by  equation  (3.15).   This  is  partially  due  to  the  fact  that  a 
different  deformation  form  was  chosen  in  reference  (30) .   However, 
a  similar  relative  relation  between  static  and  dynamic  curves, 
as  indicated  in  Figure  14  was  also  seen  in  reference  (30) .   From 
the  resulting  curves  presented  in  (30) ,  by  utilizing  the  buckling 
criterion  in  Chapter  I,  the  critical  dynamic  pressure  can  be 
obtained  with  negligible  errors  as  compared  with  the  solution. 
From  this  point  of  view,  the  result  of  the  analysis  in  (30)  may 
be  used  as  another  justification  of  the  proposed  criteriono 

3.  A  Discussion  of  the  Results 

In  both  Chapters  I  and  II,  a  qualitative  method  has  been 
used  to  discuss  the  motion  of  the  system  and  to  determine  the 
critical  condition  of  the  dynamically  loaded  structure.   It  is 
also  evident  in  these  two  chapters  that  good  agreement  exists 
between  the  result  obtained  in  this  way  and  the  result  by  the  use 
of  the  buckling  criterion  outlined  in  the  first  chapter.   This  gives, 
at  least  qualitatively,  a  justification  of  the  proposed  criterion. 
From  the  point  of  view  of  the  applications,  the  proposed  criterion 
is  subjected  to  no  restrictions  of  any  sort,  while  the  topological 
discussion  would  meet  certain  difficulties  when  the  system  is 
nonautonomous  or  of  higher  degrees  of  freedom.   However,  it  is 


55 


rather  convincing  that  the  phase  plane  method  is  suitable  for  use 
in  discussing  the  dynamic  buckling  problems.   The  accuracy  of  this 
method  may  be  restricted  by  the  fact  of  using  a  single  degree  of 
freedom  system.   The  general  behavior  of  the  motion,  however,  is 
much  clearer  as  being  plotted  in  trajectories.   Furthermore,  through 
the  examples  given  in  the  last  two  chapters,  one  can  see  the  direct 
connection  between  the  energy  method  and  the  method  of  response 
curves,  which  were  usually  employed  in  solving  the  dynamic  buckling 
problems.   This  fact  was  clearly  shown  in  Figures  7,  8,  and  9.   The 
motion  on  the  separatrix,  which  passes  a  saddle  point  corresponding 
to  a  certain  extreme  of  the  energy  level  is  the  threshold  of  the 
substantial  change  in  the  amplitude  of  the  dynamic  response. 

In  the  application  of  the  proposed  criterion  to  dynamic 
buckling  problems  considering  higher  axi symmetrical  modes  or 
unsymcietrical  forms  of  deformations,  a  suitable  characteristic 
deformation  has  to  be  chosen.   One  of  the  examples  is  to  take  the 
mean  deflection  of  the  shall  as  the  authors  of  reference  (8)    did, 

A  more  general  problem  is  that  of  rectangular  loadings 
characterized  by  the  application  of  pressure  q  at  time  t  =  0, 
which  is  held  constant  for  a  time  duration  ^  t  and  then  suddenly 
removed.   This  would  require  the  simultaneous  solution  of  two 
nonlinear  differential  equations  of  the  following  form: 
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0  <  T  <  T  . 


7j(o)    =  7j(o)     =  O     ', 


(3.26) 


?:^  "Z:  <  '^   >  (3.27) 

where  c  is  the  dimensionless  form  of  the  quantity  A  t,  according 
to  the  definition  in  (2.4),   The  stability  problem  of  this  system 
may  not  be  solved  without  having  a  general  solution  of  (3o26)o 
It  should  be  noted  that  it  is  not  possible  to  give  a  general 
e3q)ression  for  the  response  of  the  load  q  in  (3.26)  in  a  nontrivial 
form.   This  point  has  also  been  mentioned  in  Section  2,  Chapter  2. 
It  is  for  this  reason  that  only  discussion  of  obtaining  the 
results  will  be  given  in  the  following.   It  is  also  understood  that 
the  shell  buckling  occurs  after  the  load  is  removed,  i.e.,  the  time 
duration  of  the  applied  loading  is  sufficiently  small.   Therefore, 
the  loss  of  stability  of  the  system  is  largely  due  to  equation 
(3.27).   The  present  problem,  by  virtue  of  the  above  discussion, 
has  the  same  characteristics  as  the  problem  that  has  been  con- 
sidered in  Chapter  II;  in  fact,  the  latter  is  merely  a  limiting 
case  of  the  former. 

Much  as  we  have  done  in  the  previous  sections,  this 
problem  also  can  be  phrased  in  the  language  of  the  topological 
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method.   After  doing  this,  the  condition  of  instability  can  be 
formulated  in  a  straightforward  manner.   The  motion  of  the  shell 
under  the  action  of  the  said  dynamic  loading  is  described  by  both 
equation  (3.26)  and  (3,27).   Specified  clearly  in  these  two  equa- 
tions, the  motion  will  be  on  the  trajectories  of  (3.26)  until 
Z     =   Z   .      Immediately  after  Z  =  T   ,    the  motion  of  the  shell  is 
on  the  trajectories  characterized  by  equation  (3.27).   It  is 
obvious  that  the  displacement  "^  at  Z.  =  Z  is  the  common  solution 
of  both  (3.26)  and  (3.27).   The  critical  condition  of  the  struc- 
ture will  correspond  to  the  following  situation  in  the  phase  plane: 
the  response  of  q  and  its  time  derivative  in  equation  (3.26)  at 
time  Z    =  Z,  i.e.,  Viz)    and  ^(  t  )  ,  vs^hich  gives  the  initial 
conditions  for  (3.27),  will  put  the  motion  of  the  structure  on  the 
separatrix  of  system  (3.27).   We  shall  discuss  this  matter  as 
follows:  when  z  -    Z.    *    the  equation  of  the  trajectory  is 

(3.28) 

-V  2 

when  Z  >  I   ,  the  equation  of  the  trajectory  of  the  motion  reads: 


2      r  2,3 

fc,^-  =   -[(276  +  3025  0^)  f     -    ?.^4  7" 

+  1  o^U  i"" ""  -   C    1 


(3.29) 


Refer  to  equation  (3.7). 

2 
Refer  to  equation  (2.9) « 
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The   continuation  property  of   the   system  at    X    =    Z    ,    requires: 

%  (z)  =  y''(z)        ,    ^  (T)  =  ^''■(r)  . 

Therefore,  we  have 


(3.30) 


by  comparing  equation  (3.28)  and  (3.29) >   The  critical  condition 
of  the  system  under  load  q  now  depends  entirely  upon  the  value  of  C. 
Reference  is  made  to  equation  (2.18),  which  defines  the  equation 
of  the  separatrix;  the  condition  of  instability  of  the  system 
characterized  by  equations  (3.26)  and  (3.27),  or  alternatively, 
equations  (3o28)  and  (3.29)  then  turns  out  to  be 


Cs 


i,  ^(r)  =         V^  '2^  9    , 


(3.31) 


where  C  has  been  defined  in  (24.18)  and  (2.19).   As  an  example, 
when  9  =  0.26,  according  to  the  previous  analysis  in  Chapter  II, 
Cg  =  0.5183o   Therefore, 

^^^  '^^(f)   =    0  h23£6     .  (3.32) 

It  should  be  noted  that  in  equation  (3.32),  "H  (z  )    is  also  a  func- 
tion of  q.   This  is  obvious  as  shown  in  equation  (3.26) „   Since 
no  analytical  form  of  the  solution  of  (3.26)  can  be  given,  further 
discussions  would  require  a  great  number  of  calculations.   A 
procedure  for  determining  the  critical  pressure,  q^r*  ^^  suggested 
as  follows: 
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A.  A  time  duration  z     was  preassigned,  based  on  the 
external  iinpulsive  loading, 

B.  Assign  also  a  series  of  numerical  values  for  the 
loading  q  in  (3.26);  these  values  are  arranged  in  an  ascending 
order  of  magnitude  and  vith  sufficiently  small  increment.   By 
taking  9  as  a  certain  value,  e.g.,  0.26,  corresponding  to  each  q, 
every  equation  in  the  form  of  (3.26)  can  be  integrated  either 
analytically  (in  terms  of  Jacobian  Elliptic  Functions),  or 
numerically.   Therefore,  the  response  of  each  q  at  time  Z    =  Z    , 
i-e-)  '^   (r)  can  be  found. 

Co  The  critical  pressure  <\^^  is  the  one  which  satisfies 
the  condition  (3.31),  or  when  Q  =  0.26,  the  condition  (3.32)  is 
satisfied.   It  is  tedious,  yet  straightforward. 


CHAPTER  IV 


BUCKLING  OF  A  CLAMPED  SHALLOW  SPHERICAL  SHELL  UNDER  A 
UNIFORMLY  DISTRIBUTED  PRESSURE  LINEARLY  INCREASED 
WITH  TIME^ 


1 .   The  Solution  of  the  Problem 

Equation  (2.1)  in  Chapter  II  can  be  written  in  the 
following  form: 


i^)   -   ^   [i2-^^  -^W^»))< 


(4.1) 


2 
If  the  following  definitions  are  employed, 


(4.2) 


Refer  to  Figure  1  for  the  geometry  of  the  shell. 

2 
Refer  to  equations  (2.2)  and  (2.4). 
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Then  (4„1)  can  be  written  as 


+  Z  log  71^  1   _     3J  n  /  i^° 


(4.3) 


If  the  load  q  is  a  function  which  is  linearly  increasing  with 
time  t,  ioB.,  the  form  as  shown  in  Part  C  of  Figure  6,  then  it 
can  be  represented  by  the  expression: 

^   =    Q  <■    ,  (4.4) 

where  Q  is  the  pressure  increasing  rate  and  has  the  dimension 
psi.   Following  reference  (_1 )  >  we  shall  use  a  "transformed"  time 
variable  -j  *, 

^••-  =//p,   =     <%'*     ■ 

where  p^  is  the  nondimensional  critical  load  for  a  complete  shell 
given  by  classic  linear  theory  and  has  the  following  form: 
2 


n  = 


V3C'-V^} 


for  ^  =  0.3,  p  =  1.21;  Q  is  defined  in  the  following: 


Q     =       ^^VeK^ 
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which  has  a  dimension:  1/sec.   By  utilizing  (A. 4)  and  (4^5), 
equation  (4.3)  reads: 


d'rj 


idz') 


7  2  r 

— ,   +  y   (2  76  +  5.02^  eM  -7  -  ^3  1  7)' 


'       16^^^^.   (4.6) 


where 


(A.  7) 


,  a  nondimensional  quantity. 


It  is  clear  that  the  character  of  the  solution  of  (4.6)  depends 
entirely  upon  the  two  parameters  ^    and  9  •   We  also  like  to 
mention  here  that  the  transformed  time  variable  x*  ^^   the 
"dynamic  overload  factor";  i  e.,  the  ratio  of  the  critical  dynamic 
load  to  the  corresponding  critical  static  load  of  a  complete 
sphere,  which  has  been  defined  in  equation  (4.5).   Therefore,  the 
response  curve  obtained  by  integration  (46)  is  actually  the  load- 
deflection  relation  of  this  problem.   Let  us  assume  the  initial 
conditions  as  follows: 

rj  (  J-  =  0  )  =  0 

(4„8) 
which  imply  an  initially  undistrubed  shell. 
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The  nonautonomous  system  (4.6)  with  initial  conditions 
(4o8)  is  best  solved  by  a  numerical  method.   Different  values  of 
the  rates  of  dynamic  loading  and  geometrical  shape  of  shells  have 
been  selected  to  substitute  the  parameters  ^'    and  9  in  (4.6),  and 
response  curves  were  obtained  by  integrating  the  equations 
numerically  on  the  University's  IBM-7090  computer. ■'•  Dynamic 
buckling  loads  were  determined  by  use  of  the  criterion  proposed 
in  Chapter  1=   The  static  load  response  curves  were  found  by 
dropping  the  inertia  term,  i.e.,  d  '7/d-j;*   from  (4.6),  and  the 
typical  form  of  those  curves  was  shown  in  Figure  15. 

In  Figure  15,  the  form  and  the  nature  of  the  response 
curves  are  very  similar  to  their  counterparts  in  reference  (I). 
For  a  rapidly  applied  load  with  larger  Q(e.g.,  ^=   0.3),  deflec- 
tion increases  slowly  at  the  beginning  and  has  a  vigorous  change 
at  the  time  of  bucklingc   For  a  certain  shell,  i  e.,  9  ,    E, 
etc.,  are  fixed,  the  faster  the  rate  of  increasing  the  dynamic 
load,  the  higher  the  dynamic  overload  pressure  will  be.   It  is 
also  seen  in  the  same  figure,  for  a  very  slow  rate  of  load 
(e.g.,  "Ijy  =  100),  the  buckling  of  the  shell  approaches  the  static 
case  as  it  should  be  and  the  "creep  pheonomenon"  strongly  indi- 
cates that  loss  of  stability  is  of  the  "classic  type."   Followed 
by  several  cycles  of  oscillation,  the  dynamic  curves  for  W=    100 
converge  to  the  static  curve. 


^A  "Runge-Kutta  method"  was  employed.   The  technique  of 
this  method  is  found  in  (22). 
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Another  set  of  curves  was  presented  in  Figure  16,  based 
upon  various  geometrical  parameterso   As  a  different  feature  from 
the  static  case,  it  is  found  that  the  critical  load  decreases 
monotonically  with  the  increasing  Q   values  (or  decreasing  in 
shallowness).   As  shown  in  Figure  18,  no  relative  minimum  corre- 
sponding to  a  certain  d   value  seems  to  exist  as  it  usually  does 
in  the  static  case.  (Refer  to  Figure  5.)  Also  as  indicated  in 
Figure  16,  the  curves  move  toward  the  left  as  the  G    value 
increases,  corresponding  to  a  decrease  in  critical  pressure. 
However,  for  a  very  shallow  shell,  i.e.,  a  sufficiently  large 
value  of  Q   ,    e.g.,  9  =  10  (  A  xs  0„806),  as  indicated  in 
Figure  16,  the  curve  does  not  follow  the  above  argument  and  falls 
to  the  right  of  the  curves  with  Q    values  smaller  than  10= 
Because  this  curve  remains  at  very  small  deflections  at  a  very 
large  pressure,  it  is  clear  that  the  failure  of  this  structure 
will  not  be  by  buckling   As  is  well  known  in  hydrostatic  loading 
problems,  buckling  will  not  occur  for  a  very  shallow  shell  which 
has  a  geometry  close  to  a  circular  plate.   This  is  also  observed 
in  this  dynamic  loading  problem.   The  limiting  value  of  ^  ,  i.e., 
the  largest  value  of  9    for  buckling  to  occur,  has  not  been  found. 
It  is  the  feeling  of  the  author  that  the  limiting  value  of  6 
depends  also  upon  the  rate  of  the  loading,  i.e.,  the  parameter  t//o 
Therefore,  the  general  answer  may  not  be  found  without  considerable 
costly  computations.   A  further  discussion  on  this  matter  will  be 
given  later. 
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Figure  19  shows  the  relation  between  the  critical  dynamic 
overload  factor  and  the  critical  central  deflection  of  the  shell, 
i.e.,  the  deflection  at  the  time  of  the  buckling.  In  this  figure 
we  observe  the  rate  of  change  of  the  upper  ^   is  decreasing  as 

J*   increases.   This  phenomenon  can  be  explained  as  due  to  the 
development  of  membrane  stresses  which  usually  play  an  important 
role  in  the  large  deflection  theory  of  plates  and  shells.   Further- 
more, these  curves  seem  to  approach  asymptotically  to  different 
limitSo   These  limits  depend  solely  upon  the  geometry  of  the 
shell;  the  shallower  the  shell  is  the  higher  the  limit  of  the 
upper  7^j,  values  will  be.   However,  there  exists  a  particular 
shape  of  the  shallow  shell,  which  corresponds  to  0  «■  0o35.   For 
any  other  shape  of  the  shell  shallower  than  this  value  (i.e., 
6  >      0,35),  the  critical  deflection  decreases.   This  is  to  say, 
there  exists  a  maximum  for  the  upper   ^^j,  -  9        curves  as  shown 
in  Figure  20.  We  do  not  know  whether  this  value  (  0  ^  0.35)  is 
the  common  maximum  point  for  all  values  of  '(^  ;  it  would  require 
much  computation  to  answer  this  question.   We  are  rather  interested 
in  the  significance  of  the  existence  of  such  a  point.   Intuitively, 
it  is  reasonable  to  believe  that  a  shallow  shell  would  permit  more 
severe  deflection,  as  compared  to  the  height  of  its  raise,  than 
a  nonshallow  one  wouldo   This  is  noted  as  being  true  for  all  cases 
where  Q   <  0.35  in  Figure  20.   However,  a  contradiction  arises  for 
6     >      0.35  (or,  A  <  4.31) o   Two  possible  hypotheses  may  be 
provided:  first,  for  all  Q  >   0.35  there  is  no  dynamic  buckling 
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occurring  because  of  the  shallowness  of  the  shell  geometry;  second, 
for  9  ~  0.35  the  deflection  of  the  shell  has  the  highest  sensi- 
tivity, or  the  weakest  shape  with  respect  to  the  dynamic  load^ 
Therefore,  one  may  think  that  this  value  of  O    may  have  something 
to  do  with  the  size  of  the  buckling  region  of  a  complete  sphere. 

Both  Figures  19  and  20  indicate  that  if  we  want  to  select 
a  certain  value  of  characteristic  deflection  as  the  measure  of  the 
buckling  of  the  shell,  then  it  may  be  necessary  to  take  a 
different  value  of  this  quantity  for  different  geometrical  shapes 
of  the  shell  as  well  as  for  various  dynamic  loads. 

Comparing  Figures  17  and  21,  we  are  able  to  conclude  that 
the  effect  of  the  load-rate  parameter  '^   on  the  critical  dynamic 
load  is  much  more  significant  than  that  due  to  the  geometrical 
parameter  Q   .   We  mentioned  that  the  critical  dynamic  overload 
factor  Zcr*  ^^  ^  function  of  two  parameters  (variables)  Y  and 
Q   ,  i.e.,  it  is  a  surface  in  the  (  "»i/,  9  ,   X  *)  space.   Based 
upon  the  curves  in  Figures  17  to  21,  we  have  the  approximate 
equations  for  the  surface:  (For,  V  =   0.3) 

[zc;]^  -    •  ( it^-° ')[i.^i  4  2.08- e''^]- 1.32 1 


Icr 


90 


(21+w     )  1^1  +  <^-^9  e       -  J.09 


(4.9) 
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We  found  that  in  the  range  of  0.5  i  ^  i   3.5  and  0.15  i     9   ^     0„35, 
results  obtained  from  equation  (4.9)  agree  well  with  that  solved 
from  equation  (4.6)  and  initial  condition  (4  8).   The  accuracy  of 
(4.9)  can  be  proved  only  by  experiment.   We  believe,  however,  that 
a  design  formula  for  different  values  of  Poisson's  ratio  can  also 
be  established  in  the  following  form 


icr  -   (^  +  v'^Cc  f  }  e-'^)     ,         ^^_^Q^ 


with  a,  b,  c,  g,  and  k  determined  by  experimental  testSo 

The  projection  of  the  trajectories  on  the  ^  ~  V     pl-ane 
are  shown  in  Figures  23  and  24.   According  to  the  criterion  in 
Chapter  I,     the  critical  state  is  bounded  by  the  two  points 
indicated  by  L  and  U,  i  e,,  the  lower  and  the  upper  bounds.   It 
is  seen  from  these  trajectories  that  the  deformation  ''^  increases 
very  slightly  after  passing  the  point  U  and  oscillates  about 
different  equilibrium  points  on  the  ''7-axis  with  an  increase  in 
amplitude. 

2o   A  Conclusive  Discussion  of  the  Problem 

The  nature  of  the  response  of  shallow  spherical  shells 
to  a  high-speed  dynamic  load  with  linearly  increasing  pressure 
intensity  were  found  and  represented  in  Figures  15  to  24.   We 
found  the  functional  relation  between  the  critical  load  and  the 
geometry  of  the  shall  has  a  characteristic  difference  from  the 
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static  case  and  this  is  shown  by  comparing  the  curves  in  Figure  17 
with  their  counterpart  in  Figure  5. 

Differing  from  the  result  obtained  by  Agamirov  and 
Vol'mir,  we  found  the  critical  dynamic  load  (or,  D.O.L.F.  ) 
depends  upon  the  ratio  of  the  speed  of  elastic  waves  in  the  shell 
material  and  the  product  of  the  shell  radius  with  the  increasing 
rate  of  the  intensity  of  the  external  pressure,  instead  of 
solely  upon  the  increasing  rate  of  the  load  as  presented  in  (l) . 
A  functional  relation  between  the  critical  dynamic  overload  pres- 
sure and  the  two  parameters  1//"  and  6    Vvas  formulated  in  equa- 
tion (4  y) ;  it  gives  the  critical  dynamic  load  from  the  given 
geometric  shape,  the  material  properties,  and  the  increasing 
rate  of  the  dynamic  loading. 

Another  suggestion  was  also  offered  by  formula  (4.9). 
Let  us  return  to  the  definition  of  "l/^  in  (4.7),  i.e., 

-^  =     ^-f  ^A5)      . 

which  is  a  dimensionless  quantity.   The  requirement  of  similitude 
is  very  conveniently  furnished  by  the  quantity  ^/  .      For  an 
example,  if  we  choose  the  same  geometric  parameter  Q    for  the 
model  and  the  prototype,  then  one  can  determine  the  nature  of 


Dynamic  overload  factor. 
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dynamic  buckling  of  the  prototype  under  a  very  high  rate  of 
dynamic  loading,  i.e„,  a  very  large  Q,  from  the  test  of  a  model 
with  lower  Young's  Modulus  and  a  relatively  lower  rate  of  the 
dynamic  load,  provided  they  have  the  same  value  of  ^, 

Discussion  of  the  buckling  criterion  was  also  made.   It 
is  seen  from  the  results  that  both  the  increasing  rate  of  the 
dynamic  load  and  the  geometrical  shape  of  the  shell  have  influ- 
ences on  the  critical  deflection  of  the  structure.   Therefore, 
the  criterion  proposed  in  Chapter  I,  which  permitted  "the  measure 
of  the  critical  state,"  to  change  along  with  different  dynamic 
loadings  and  geometric  shapes  of  the  shall,  has  definite  advantages 


CHAPTER  V 


BUCKLING  OF  A  COMPLETE  SPHERICAL  SHELL  UNDER  SUDDENLY 
APPLIED  UNIFORMLY  DISTRIBUTED  DYNAMIC  LOADINGS 


Introduction 

The  nonlinear  problem  of  a  complete  spherical  shell  has 
been  discussed  by  numerous  people  since  the  article  by  von 
Karman  and  Tsien  first  appeared  in  1941  (17) ,   The  progress  toward 
a  substantiated  explanation  of  the  discrepancy  between  the  classic 
linear  theory  of  Zolley  and  experiments  has  been  rather  slow. 
This  is  partially  due  to  the  lack  of  reliable  experimental  data 
on  the  buckling  load  for  a  complete  spherical  shell.   For  many 
years  Tsien' s  (31)  energy  criterion  and  the  "lower  buckling 
load"  (32)  have  been  used  to  determine  the  load-carrying  capacity 
of  a  thin  elastic  spherical  shell.   In  the  year  1962,  two  inde- 
pendent experimental  tests  contributing  to  the  problem  were 
reported  by  Thompson  (33)  and  Krenzke  (19) «   Both  of  these  two 
experiments  show  that  the  critical  load  of  a  complete  spherical 
shell  can  be  much  higher  than  the  lower  buckling  load  offered  by 
Tsien;  it  ranges  from  45  per  cent  to  70  per  cent  of  the  critical 
pressure  predicted  by  classic  linear  theory.   A  fairly  precise 
theoretical  analysis  based  upon  a  large-deflection  strain  energy 
theory  was  also  given  by  Thompson  in  reference  (33) .   Thompson 
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ended  his  research  with  the  conclusions: 

A.  The  initial  buckling  (pre-buckling)  was  seen  to  be 
classical  in  nature,  i.e.,  the  load  deflection  process  is  con- 
tinuous, 

B.  The  buckling  load  was  about  75  per  cent  of  Zolley's 
"upper  buckling  load." 

Co  The  stable  post-buckling  states  are  observed  to  be 
rotationally  symmetrical o 

Both  Thompson  and  Krenzke's  tests  suggested  further 
theoretical  investigations.   It  seems  that  more  precise  work  on 
the  numerical  investigation  of  the  governing  differential  equa- 
tions as  has  been  done  for  shallow  shells  is  very  much  to  be 
desiredo^ 

Among  all  the  works  on  the  stability  problem  of  complete 
spherical  shells,  Vol'mir's  first  approximation  (35)  remains 
unknown  to  most  of  the  western  scholars,  mainly  because  it  was  not 
a  successful  one.   He  attempted  to  investigate  the  post-buckling 
behavior  of  a  complete  spherical  shell  by  using  a  very  simple  form 
of  deformation,  i.e.,  the  form  in  equation  (2.0)  and  integrated 
the  equation  approximately  by  a  Galerkin  method.   His  assumption 
on  the  buckling  process  was  the  same  as  the  one  assumed  by  von 


^E.g.,  (16),  (7),  (27),  (38),  (39).  and  (40) 
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Karman  and  Tsien,  i.e.,  the  sphere  is  contracted  to  a  smaller 
sphere  and  forms  a  single  dimple  after  buckling   By  virtue  of 
this  buckling  mechanism,  Vol'mir  allowed  a  uniform  membrane 
stress:  ff' =  -  "  "/^A.  to  be  distributed  in  the  buckling  region  before 
the  loss  of  stability  and  considered  the  buckling  region  as  a 
shallow  spherical  cap.   In  othir  words,  the  differences  between 
Vol'mir' s  method  and  Karman-Tsien' s  is  that  Vol  mir  used  a 
Galcrkin's  method  based  upon  a  variation  of  the  governing 
differential  equations  and  also  considered  the  variation  only  of 
a  single  parameter,  which  is  the  central  deflection,  while  Karraan- 
Tsten  based  their  method  upon  the  variation  of  the  total  energy 
with  two  parameters.   In  a  limiting  case  of  small  deflection 
theory,  Vol'mir  found  a  buckling  load  twice  as  large  in  magnitude 
as  predicted  by  the  classic  linear  theory  and  failed  to  obtain 
a  lower  buckling  load. 

In  Vol 'mir' s  approximation  method,  we  find  certain  ques- 
tions which  may  have  led  to  his  failure  in  obtaining  an  approx- 
imate solution.   First,  he  considered  that  buckling  occurs  at  the 
transition  of  the  membrane  state  and  the  bending  state  of  the 
shell;  therefore,  his  effort  was  concentrated  on  finding  a  post- 
buckling  load.   Second,  the  assumption  of  a  uniform  membrane 
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stress  in  the  buckling  region  before  the  loss  of  stability 
requires  an  abrupt  change  of  stress  state,  which  seems 
impossible  under  a  continuous  loading  process.   If  the 
structure  is  perfect  in  geometrical  shape,  this  change  of 
stress  distribution  would  require  different  equations  to 
describe  the  equilibrium  of  the  shell  and  result  in 
mathematical  difficulties  of  integrating  the  differential 
equations.   It  seems  to  us  that  Vol'mir's  method  may 
be  restudied  using  a  different  contemplation  of  the  buckling 
mechanism  of  a  complete  spherical  shell. 


1 _   On  a  New  Mechanism  of  the  Buckling  of  a  Complete  Spherical 
Shell 


We  shall  assume  that  the  buckling  of  the  shell  follows 
a  possible  mechanism  which  permits  the  transition  from  a  membrane 
state  of  stress  to  a  bending  stress  state  occurring  in  a  con- 
tinuous fashion  and  the  transition  occurs  before  the  loss  of 
stability  of  the  shell.   This  is  described  in  the  following 
paragraphs  c. 

A.   When  the  external  pressure  q  is  much  less  than  the 
critical  value,  the  shell  contracts  to  a  slightly  smaller  sphere; 
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as  shown  in  Figure  25,  Part  A,  the  original  shell  contracts  to 
a  sphere  with  radius  equal  to  Rq  -  ^  Ro   As  q  increases,  AR 
increases  and  a  significant  change  in  curvature  occurs  because 

of  the  change  in  radius. 

B.  Let  us  take  another  assumption  that  the  structure  has 
a  resistant  nature  against  the  higher  pressure  and  has  a  tendency 
to  resume  its  original  curvature.  Based  on  experimental  evidence 
that  spheres  form  a  single  dent  after  buckling,  we  thii\k  that  the 
resumption  of  the  curvature  starts  at  a  small  region,  or  we  might 
say  at  a  point.  The  effect  of  the  resumption  of  the  curvature 
from  a  larger  one  to  a  smaller  one  has  introduced  a  pure  moment, 

which  will  be  in  the  same  direction  as  that  caused  by  the  external 
pressure.    In  other  words,  we  consider  that  the  initiation  of 
bending  stress  in  the  shell  is  due  to  the  imperfect  nature  of  the 
structure;  however,  the  external  pressure  will  certainly  help  to 
increase  the  magnitude  of  the  bending  stress  and  build  up  the 
inward  deflections » 


Refer  to  Figure  25,  Part  C. 
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Co   The  existence  of  such  a  single  point  for  the  first 
resun^Jtion  of  the  original  curvature  may  be  explained  as  due  to 
the  "imperfections  "   Let  us  assume  a  spherical  shell  with  perfect 
geometrical  shape  all  around  except  a  very  small  hole  at  point  0. 

(The  advantage  of  the  assumption  of  a  small  hole  is  that  we 
do  not  have  to  make  any  other  assumptions  on  the  form  of  the 
imperfections  a)   When  the  original  shell  contracts  to  a  smaller 
one  (refer  to  Figure  25)  so  that  A  moves  to  A',  and  0  to  0', 
the  small  hole  is  contracted  to  an  infinitesimal  one.   Under 
such  a  condition,  the  membrane  stress  at  0'  is  certainly 
zerOo 

Do   If  we  allow  the  existence  of  such  an  infinitesimal 
hole  at  the  point  0  in  Figure  25,  then  the  bending  state  is 
inherent  in  the  problem  itself.   As  shown  in  Figure  25,  Part  C, 
in  the  immediate  vicinity  of  0',  the  situation  is  very  similar 
to  a  clamped  circular  plate  with  a  central  hole.   The  idea  of 
"boundary  layer"  may  be  best  fitted  into  this  particular 
circular  region;  outside  of  this  region,  a  pure  membrane  state 
remains.   When  the  external  pressure  q  increases,  this  circular 
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region  dilates  (or  the  thickness  of  the  boundary  layer  increases) 
and  forms  the  buckling  region  after  the  loss  of  stability  of 
the  shell „ 

By  virtue  of  the  a  bove  described  buckling  mechanism, 
we  arrive  at  the  conclusion  of  the  existence  of  a  boundary 
layer  at  the  vicinity  of  a  point  0' .      In  this  region,  both 
membrane  and  bending  stresses  exist  at  the  time  of  the  stabil- 
ity of  the  shell.   We  are  interested  in  the  distribution  of  the 
membrane  stress  in  the  boundary  layer  region  during  the  load- 
ing process.   As  we  have  mentioned,  the  stress  at  0'  is  zero 
and  outside  of  this  region  the  shell  maintains  a  raomentless 
state  with  a  membrane  stress   OV  =  -  o"o /2h  .      By  referring 
to  the  stress  distribution  in  a  bent  clamped  circular  plate, ^ 

a  reasonable  assunqjtion  in  the  boundary  layer  region  will  be 

2 
a  parabolic  variation,   i.e., 


^Refer  to  Q4)  ,  pp.  54-63 » 
^Refer  to  Figure  25,  Part  Bo 
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We  shall  analyze  a  nonlinear  problem  of  the  loss  of 
stability  of  a  complete  sphere  by  taking  the  buckling  region  as 
a  clamped  shallow  spherical  segment  with  a  nonuniform  membrane 
stress  in  the  form  of  (5.1)  distributed  in  the  middle  surface 
before  its  loss  of  stability. 

2„   Buckling  of  the  Shell  under  a  Static  Load 

We  shall  take  the  buckling  region  of  a  complete  spherical 
shell  as  a  shallow  spherical  segment  clamped  along  a  circular 
boundary,  as  shown  in  Figure  1,  Part  B.   From  the  discussion  in 
Appendix  II  for  shallow  spherical  shells,  we  have  the  governing 
differential  equations  for  such  a  shell  under  a  uniformly 
distributed  static  load  q  (p„s  i.)  in  the  following  form: 

The  equation  of  equilibrium, 

(5.2) 
and  the  compatibility  equation, 


(5.3) 


where 


^2  =      d_l  LJ_ 


Refer  to  equations  (Ao29)  and  (A. 30)  in  Appendix  II. 
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Since  we  have  restricted  ourselves  in  the  problem  of  axi- 
synmetrical  deformations,  then  equation  (2.0)  can  be  used  again 
as  a  first  approximation  of  the  deflections  in  the  buckling 
region;  thus,  we  have 

2 


W  = 


«['-('A/)]     .  (5.4) 


Substituting  (5o4)  into  (5.3)  and  integrating,  we  have 


d 


_  e^  1^0  r  2(2-;))/  r 


6R 
which  is  the  condition  that  the  strains  or  stresses  in  the  middle 

surface  due  to  large  deflections  have  to  satisfy.   However,  it 
should  be  remembered  that  there  is  a  membrane  stress  already  in 
the  middle  surface  due  to  the  contraction  effect  of  the  rest  of 
the  spherical  shell  outside  of  the  boundary  layer  regiono   By 
using  the  relation  between  the  stress  function  <p    and  radial 
stress  (Ty    ,  and  after  including  equation  (5.1),  the  above  equa- 
tion of  the  compatibility  condition  of  deformations  in  the  middle 
surface  of  the  shell  reads 

r  .1 


6flo  L  (>~^^  ^""o^        ^^'^    ^"'^  JJ    2h    Ui 

(5.5) 
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Substituting  (5.5)  into  (5.2),  the  equation  of  equilibrium  yields 
the  following  form 


Gr 


Ta 


6r< 


£^  ro  r  2(2-1^)/  r 


-(.^)^ 


Mfj^ 


2B 


2h  (^'o^ 


i'f'-^jj 


(5.6) 


We  shall  use  Galerkin's  method.  This  method  allov.s 
equation  (5.4)  and  (5  5)  to  be  the  approximate  solutions  of 
equations  (5o2)  and  (5.3),  provided 


JJ  A 


(5.7) 


is  satisfied.   In  (5.7)  "A"  represents  the  area  of  the  circular 

1 
region  with  radius  r  .    It  should  be  noted  that  r   is  not  a 
o  o 

constant;  it  is  the  thickness  of  the  boundary  layer,  which 
varies  with  the  external  pressure. 

After  performing  the  above  int'  gration,  we  find  that  the 
central  deflection  S  has  to  satisfy  the  following  equation, 
which  describes  the  equilibrium  conditions « 


■Refer  to  Figure  25. 
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[(.->'")  C      ^     K   1  24  (.-J       12  \  j     ^ 


RoTo'  I 


^n  -275  >< 

k9  (.i-v) 


(,219 

nou. 


]( 


1  r  ^-^ 


30 

/3 


r  = 


(5.8) 


For  ^=  0.3,  equation  (5.8)  reads 


12/  ^^   ,   276 


)^ 


lo_62  2 


(5.9) 


As  in  the  case  of  shallow  shells,  we  shall  employ  the  approx- 
imations and  diraensionless  quantities  described  in  equations 
(2.2)  and  (2,4),  i.e.. 


(5.10) 


It  should  be  noted  that  among  these  quantities,  z   is  neither  a 

o 

fixed  constant  nor  a  given  value  as  in  the  case  of  shallow  shells. 
It  depends  upon  the  size  of  the  buckling  region  as  indicated  in 
(5.10).   For  the  same  reason,  Q   should  be  treated  as  a  parameter 
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in  the  sequel,  and  it  will  be  used  as  a  measure  of  the  size  of  the 
buckling  region,  or  the  thickness  of  the  boundary  layer,  to  replace 
the  variable  r^. 

By  utilizing  (5.10),  equation  (5.9)  takes  the  following 
form 

f    =    |[%   +  ^  e?]] -[(276 +  3.025  9^)'rj    -   5  31   ^^ 

+  2108-  rj^J        =  0     .       (5«11) 

this  is  the  load-deflection  relation  when  the  shell  is  under 
static  equilibrium^ 

Let  us  take  this  opportunity  to  discuss  the  nature  of 
the  load-deflection  relation  and  its  associated  stability  proper- 
ties.  A  typical  curve  of  equation  (5.11)  is  in  the  form  of 
curve  1  as  shown  in  the  figure  on  the  following  page.   Curve  1 
constitutes  three  branches:  the  unbuckled  stable  branch  OA,  the 
unstable  branch  AB,  and  the  buckled  stable  branch  BA' .   Instead 
of  calling  point  A  and  B  the  bifurcation  points  or  branch  points, 
we  shall  directly  call  them  the  critical  points.   The  feature  of 
the  loss  of  stability  is  such  that,  during  the  loading  process, 
the  equilibrium  position  of  the  structure  moves  from  0  to  B'  and 
then  to  A  in  a  continuous  and  monotonously  increasing  fashion; 
any  slight  increment  of  the  pressure  at  the  equilibrium  position 
at  A  would  cause  a  sudden  and  large  increase  in  deflection,  which 
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brings  the  equilibrium  position  from  A  to  some  point  above  A' 
on  the  buckled  stable  branch   Therefore,  the  equilibrium  condi- 
tion at  A  is  certainly  a  "critical"  situation,  and  the  load 
corresponding  to  the  equilibrium  condition  at  A  on  curve  1  deserves 
the  name  of  "critical  load."  In  the  static  loading  analysis,  we 
shall  permit  only  one  type  of  buckling  of  the  shell  that  is  due 
to  the  loss  of  stability  after  passing  the  point  A;  therefore, 
the  critical  load  at  point  A  is  also  the  buckling  load  of  the 
structure.   The  significance  of  the  equilibrium  situation  at  the 
critical  point  B  defines  the  equilibrium  condition  where  the 
"outward  snap"  of  the  shell  occurs.   For  the  purpose  of  emphasis, 
we  shall  repeat  the  argument  that  equilibrium  condition  at 
point  B  hini   no  significance  to  the  instability  of  the  shells,  if 
a  classic  buckling  criterion  is  used.   We  would  like  also  to 
point  out  that  it  is  incorrect  when  we  have  a  P-  d  curve  in  the 
form  of  curve  1;  this  then  implies  the  necessity  of  using  an 
"energy  criterion  "   The  main  difference  between  these  two 
criterions  is  the  method  of  determining  the  buckling  load.   The 
classic  criterion  defines  the  buckling  load  by  having  a  hori- 
zontal tangent  at  the  critical  point,  i.e.,  the  buckling  of  the 


We  use  the  definition  given  by  Kaplan  and  Fung  in 
reference  (16) „ 
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shell  is  solely  due  to  the  loss  of  stability  and  the  load-deflection 
curve  usually  is  in  the  form  of  curve  II  in  the  figure  referred  to. 
On  the  other  hand,  energy  criterion  permits  a  "jump"  from  the  state 
at  B'  to  B  during  the  inward  deformation  process;  thus,  a  lower 
buckling  load  corresponding  to  the  pressure  at  the  equilibrium 
state  at  B  is  defined.   Curves  I  and  II  indeed  represent  two 
different  types  of  instability.   According  to  Biezeno  and  Graramel  (b) , 
the  instability  represented  by  curve  I  is  called  "transitional 
instability"  and  the  other  is  called  "complete  instability,"  which 
does  not  have  the  monotonously  increasing  branch  beyond  the  range  of 
instability.   An  example  based  upon  a  classic  criterion  was  given  in 
reference  (6),      where  the  load  deflection  relation  was  in  the  form 
of  curve  I.   We  shall  investigate  the  instability  of  the  shell,  basing 
the  investigation  upon  a  classic  criterion.   The  condition  of  equili- 
brium corresponding  to  point  A  in  the  above-mentioned  figure  will  be 
referred  to  as  the  "unstable  equilibrium  position"  or  the  "critical 
position"  because  it  defines  the  loss  of  stability  of  the  system  as 
well  as  the  buckling  of  the  shell. 

Let  us  return  to  equation  (5.11).  Geometrically,  it  repre- 
sents a  one  parameter  family  of  curves  in  the  ^ -q  plane;  6  is  the 
parameter.   For  each  8  value,  equation  (5.11)  shows  a  possible 


Refer  to  reference  (6),  pp.  484-496. 
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load-deflection  relation  during  the  buckling  process.   Let  us  say 

that  the  true  P-  S   relation  during  the  buckling  of  the  shell  will 

be  the  one  with  0=  ^cr'  ^^^  ^^  ^   form  similar  to  curve  I,  which  has 

been  discussed  above.   We  shall  define  the  6   in  the  following 

cr  ^ 

fashion:   0   will  make  the  system  reach  its  "unstable  equilibrium 
position"  with  the  smallest  value  of  q„   By  the  fact  that  the  state 
of  unstable  equilibrium  corresponds  to  a  relative  maximum  position 
on  the  q( 7 )  curve  mathematically,  the  problem  of  finding  Q       is 
equivalent  to  seeking  a  least  maximum  for  the  family  of  curves  in 
equation  (5.11). 

After  calculations,  the  locus  of  the  "unstable  equilibrium 
position"  was  found  to  be  as  shown  in  Figure  26.   It  has  a  relative 
minimum  when  6=   0.548.   Therefore,  we  have 

0cr  -  0-5^8"  (5^2) 

As  9  =   0   =  0.548,  equation  (5.11)  takes  the  following  form: 

0       _         3  66g-7f  rj    -  ^31  rj^    -¥   2lo^  ^^ 

0-^Ug  (  I  31^    -h    olS   r])  .  (5.13) 

This  equation  describes  the  curve  shown  in  Figure  26,  which  repre- 
sents the  "best  possible"  relation  between  the  load  and  deflection 
in  the  sense  of  yielding  a  smallest  critical  load.   It  posses  a 
maximum  at  ^  =  0.418  corresponding  to  the  critical  position  and  a 
minimum  at  f]  =   1.22.   Corresponding  to  7 ~  0.418,  we  have  the 
critical  load  from  equation  (5.13): 
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i        -  0.&2I  (5.1A) 

0  cr  1 

which  is  about  68  per  cent  of  Zolley's  result  based  on  a  classic 
linear  theory  and  matches  very  well  with  both  Thompson  and  Krenzke's 
experimental  results.   It  also  should  be  noted  that  corresponding 
to  the  minimum  position  on  the  ^ -q  curve  there  is  a 

%        =    0-32    ,  (5,15) 

()  oa 

which  is  about  26  per  cent  of  the  result  of  the  classic  linear  theory. 
As  we  have  mentioned,  this  load  corresponds  to  the  outward  snapping 
load  in  the  unloading  process.   The  experimental  test  by  Thompson 
gave  an  outward  snapping  load  for  shells  with  (R^/h)  ~  20  of  the 
magnitude  about  22  per  cent  of  the  linear  classic  result.   The  change 
of  volume  during  the  loading  process  can  also  be  obtained  from  the 
analysis  and  has  the  following  form: 

Av    =    211  (i-v)  —     +     yrj  ^     , 

The  first  term  in  the  above  equation  was  due  to  the  membrane  contrac- 
tion, while  the  second  term  was  the  volume  developed  by  inward 
deformation  in  the  buckling  region  and  was  found  by  the  following 
equation: 


Krenzke  found  experimentally,  q^j.  =  0.84;  refer  to 
reference  (19) . 
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^0  r2n 
Let  us  define  a  dimensionless  change  of  volume  in  the  form: 

then,  from  the  above  equations,  we  have 

where  Q       =   0.548,  defined  in  (6.12),  and  q,  f]   were  also  given 
previously.   We  shall  use  the  subscript  "cl"  to  indicate  quantities 
corresponding  to  the  result  obtained  from  the  classic  linear 
theory,  i „e. , 

q    _     2 


,2 


Therefore,  for  V  =  0.3,(R/h)  =  20,  we  obtain  the  following  expression: 

Together  with  equation  (5.13),  the  relation  between  (  ^^  /ji\/   -,)    and 
(q/q^]^)  can  be  found;  it  is  given  in  Figure  27.   As  far  as  the 
critical  loads  are  concerned,  the  theoretical  result  is  qualitatively 
good  as  compared  with  experiment.   The  rate  of  increasing  of  pressure 
in  the  post-buckling  region  was  seen  to  be  faster  than  the  experi- 
mental results  given  by  (33) « 


3.   Buckling  of  the  Shell  under  a  Suddenly  Applied  Constant  Pressure 


Let  us  assume  that  Reissner' s  simplified  theory  on  transverse 
vibration  of  a  thin  shallow  elastic  shell,  i.e.,  the  inertial  forces 
in  the  middle  surface  are  neglected  as  being  small  compared  to  the 
transverse  inertial  force  does  also  hold  true  in  the  case  of  a 
complete  spherical  shell.   We,  therefore,  obtain  the  equation  of 
motion  for  such  a  shell  by  adding  one  terra: 

which  is  due  to  the  transverse  inertial  effects,  to  the  right-hand 
side  of  the  equation  of  equilibrium  (5.2) «   Equations  of  compati- 
bility are  kept  in  the  same  form  as  equation  (5.3) o   By  taking  the 
same  form  of  equation  (5.4)  for  the  axi symmetrical  dynamic  deforma- 
tions of  the  shell  and  considering  the  central  deflection  varying 
along  with  time,  after  performing  a  similar  integration  of  the 
Galerkin's  functional   in  the  form  of  (5.7),  we  obtain  the  counter- 
part of  equation  (5.9),  i.e.,  the  dynamic  equation  of  a  metallic 
(  'L'  =  0.3)  complete  spherical  shell  in  the  following  form: 


d'^   ^   ^  F  /  '2  '  ^^  ^  2.7(,   ,  ,      10. 6Z         ,2 


(5ol7) 


C  M  -   Ph  t  ^   ^  r,^  ^  J. 


1 

Refer  to  Figure  6,  Part  B. 
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If  the  following  dlraensionless  quantities,  which  have  been  employed 
in  the  case  of  a  shallow  spherical  shell  are  used-'^ 


^-  i     '     ^=  h  '     ^  -  I 


^-  ^fijf  •    ^'=1?. 


o 

2, 


(5.18) 
then  equation  (5ol7)  takes  the  following  nondimensional  form: 

r]    +  j(276  +  302^  e^)->7  -  ^31  rf-  -{■  2  io^   ^^  J 


(5.19) 


If  a  shell  is  under  such  condition  that  it  is  initially 
undisturbed,  then  at  time  X  =   X  q  =  0,  we  have  the  following 
initial  conditions 


71  (o)      =         0  and    7J  (O)   =    0     . 

(5.20) 

2 
Because  of  the  ^   term  appearing  in  the  nonlinear  part, 

equation  (5.19)  does  not  belong  to  any  well-known  class  of  equations 


Refer  to  equation  (2.4). 
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whose  behavior  has  been  systematically  discussedc   In  the  case  of  a 
given  load  q,  the  dynamic  response  can  be  obtained  by  integrating 
(5.19)  and  the  solution,  in  general,  in  terms  of  elliptical  functions. 

Now,  in  this  problem,  q  is  taken  as  an  unknown  parameter. 
We  are  looking  for  the  critical  value  of  q  which  will  result  in  the 
motion  of  the  system  being  unstable.   Since  the  system  in  (5ol9) 
appears  to  be  autonomous,  then  a  qualitative  discussion  of  the 
motion  is  possible  by  use  of  a  topological  method. 

Equation  (5.19)  can  be  integrated  once  when  q  is  taken  as  a 
constant,  and  the  following  result  is  obtained: 


+  3^^  ^^  -  /^^^  ^^     r  >o 


(5.21) 


When  q,  9    and  the  arbitrary  constant  A  are  assigned  certain  values, 
equation  (5.21)  is  the  equation  of  the  trajectories  in  the  phase 
plane  of  the  system  (5.19).   This  equation,  in  general,  posses 
three  (3)  singular  points  at  the  7  values  solved  from  the  following 
equations: 

^  =    !>    =    0 
i>   =     V    =    0 

(5.22) 
or  on  the  fj   axis  where 
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^  ^9  -  (216   +  J025  e'-^fe)^   +  ^3)  rj' 
-  2  109-    '7^   =   0   . 


(5.23) 


Let  ''?,   '^^,  and   ^  ,  arranged  in  an  order  of  increasing  raagni- 
tude,  be  the  roots  solved  from  equation  (5.23)  when  the  values  of 
q  and  6  are  given.   Tested  by  the  stability  theorem  (9^),   we 
shall  see  there  are,  in  general,  two  centers  and  one  saddle  point. 
Because  of  the  initial  conditions  in  (5.20),  the  trajectory  which 

passes  through  the  origin  (0,0)  of  the  phase  plane  does  represent 

2 
the  oscillatory  motion  of  the  shell.    Therefore,  the  problem  of 

finding  the  critical  q    for  equation  (5.19)  is  equivalent  to 

finding  the  q  value  which  will  make  the  separatrix  pass  through 

3 
the  origin  (0,0)  of  the  phase  plane.    Equation  (5.21)  suggests 

that  for  the  trajectory  passing  through  the  origin,  the  value  of  A 

in  (5.21)  should  be  taken  equal  to  zero.   Thus, 


Refer  to  reference  {9),    p.  317. 

2 
Refer  to  the  discussions  in  Chapter  III. 


3 
Refer  to  Figure  28. 
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^'  -    ^[7  fe  -  (Z76  +  502^6^-  I  J  0)^7 
+     35-4  ^^  -  /c5-4  -7^  J   , 

(5.24) 

Let  us  assume  that  among  the  three  roots  of  (5.23), 
corresponding  to  '^  =     ^2'    the  point:   (  "^  ,  0)  is  the  saddle 
point;  then  (5.24)  will  be  the  equation  of  the  separatrix  provided 
the  following  condition  is  satisfied: 

-%[!  oik  T]l  -  3  5/,  7/  +  (2  76  +  3 »«  e'- 1  f  e)  17^ 

|e]    =    0 


l\ 


In  case  that   7  f   0,  we  have 

-[1  osu  ^1  -  ziii  if    +  (276  +  3o25e^-  I  fe)  \ 


fsj 


1':::    =     0 


(5.25) 


It  should  be  remembered  that  7~   ^o  is  one  of  the  roots  of  (5.23) 
Furthermore,  for  the  point  (  ^o'  ^^    "-^  ^^   ^"^  unstable  critical 
point,  the  following  condition  obtained  from  the  characteristic 
equation  of  the  stability  theory  must  be  satisfied,  i.e.. 


Refer  to  Appendix  I  for  discussions. 
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[  10  62    rj^   ~  6  52U    r^^  -  ( 2lG  +   3.o25  6^  -    V4    J  6 )]   >   0, 

(5»26) 

The  problem  of  finding  the  q   has  been  reduced  to 
seeking  a  common  solution  between  the  two  algebraic  equations 
(5.23)  and  (5.25)  and  satisfies  condition  (5.26). 

Now  let  us  return  to  equation  (5.19).   There  still  remains 

the  main  problem  of  choosing  the  value  of  Q    which  characterizes  the 

very  nature  of  the  instability  as  well  as  the  magnitude  of  the 

critical  load.   In  this  section,  we  shall  present  two  methods  of 

determining  the  buckling  region  parameter  9     .   Detailed  discus- 

cr 

sions,  as  well  as  numerical  examples,  will  be  given  in  both 
methods. 

Ao  A  Static  Approach  to  the  Determination  of  9cr° 


We  shaUL  go  back  to  the  last  problem,  i.e.,  to  find  a  common 
root  between  (5.23)  and  (5.25)    After  the  following  discussion,  we 
shall  see  that  the  6  value  is  bounded;  i.ec,  it  can  assume  any 
value  but  only  in  a  certain  region. 

From  equation  (5.23)  and  (5.25),  the  difference  between 
these  two  equations  can  be  found  and  has  the  form: 


111    rj^   -  I  o^k   rj^ 


y  375  9 

(5.27) 


94 


There  exists  at  leas"t  one  common  root  between  (5.23)  and  (5.25)  if 
the  relation  (5.27)  is  eraployedo  Substituting  (5.27)  into  (5.23) , 
equation  (5.23)  reads: 


+  (2i6  +  3  02^  e^ )1  V      —       0 

J     '  ■  (5.28) 


Real  nontrivial  roots  of  (5  28)  are  solved  from  the  equation  which 
is  obtained  by  setting  the  expression  in  the  bracket  of  (5.28) 
equal  to  zero.   Solutions  of  this  equation  can  be  discussed  graphi- 
cally.  Let  us  designate 

f^  =   (  2  76   -^  3  02^   e^j.  (329) 

Solutions  of  (5.28)  are  simply  the  ^  coordinates  of  the  intersecting 
points  of  curves  f,  (V)    ^^'^   ^o  ^1^'        ^°^   physical  reasons  we  are 
interested  only  in  positive  values.   Equation  (5.29)  shows  that 
f,  ( "y )  is  a  fixed  curve  in  the  plane,  but  f   ('7)  is  a  horizontal 
line  parallel  to  the  ''^    axis,  and  its  position  depends  upon  the 
parameter  6   .   The  value  of  1/  9   measures  the  size  of  the  buckling 
region  or  the  thickness  of  the  boundary  layer  that  may  appear  in  the 


See  Figure  29. 
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form  of  a  small  dent  after  the  shell  loses  its  stability.   For  this 
reason,  the  physical  problem  is  meaningful  only  when  6    has  some 
finite  value  greater  than  zero.   Referring  to  Figure  29,  it  is 
obvious  that  there  will  be  no  real  positive  roots  of  equation  (5  28) 
other  than  the  trivial  solution  if  f„  >  4  3665.   Also  because  of 
the  lower  bound  of  the  9  value,  it  makes  the  value  of  f   bounded 
in  the  range  2.76  ^  f^  <  4.3665,  or  equivalently, 

0    <       9       <       072B5   .  (5.30) 


It  is  also  obvious,  from  Figure  29,  that  in  general  (5.28)  has  two 
real  positive  roots.   Substituting  either  or  both  of  these  two  roots 
together  with  the  corresponding  q  values  solved  from  (5  27)  into 
equations  (5.23)  and  (5.25)  will  make  these  two  equations  have  common 
solutions.   It  will  be  seen  that  only  one  pair  of  results  will  cha- 
racterize the  loss  of  stability  of  the  nature  discussed  previously; 
i.e.,  the  motion  of  the  system  is  of  the  nature  described  in 
Figure  28,   From  Figure  28,  the  trajectories  numbered  in  the  figures 
represent  the  motion  of  the  system  (5.19).   When  q  <  q   ,  the 
motion  is  stable  and  oscillates  about  a  stable  singular  point 
(  ^p  0)  with  amplitude  '^    -      When  q  =  q   ,  the  motion  is  on  the 
separatrix,  which  passes  through  the  unstable  singular  point  (  V-,  0  ) 


From  (5  27),  it  is  apparent  that  when  '7=0,  there  is 
implied  an  undisturbed  condition^ 
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As  q  >  q   ,  the  motion  is  unstable  and  a  severe  change  of  maximum 

deflection,  as  indicated  by  comparing  ^  and  V   ,    vill  result  in 

p       m 

the  buckling  of  the  shell.   The  second  pair  of  the  possible  results, 
as  tested  by  (5  26) ,  is  stable.   However,  there  exists  the  possi- 
bility that  the  oscillatory  motion  about  one  stable  singular  point 
may  jump  into  another  state  of  oscillatory  motion  about  a  different 
stable  singular  point.    Because  of  the  sudden  and  severe  change  of 
maximum  deflection,  buckling  of  the  structure  may  also  occur,  and 
it  may  be  more  important  than  the  first  type  in  engineering  appli- 
cations.  The  method  of  finding  the  q   's  based  upon  this  analysis 
is  illustrated  in  the  folloving  example. 
When   0  =  0.  44  ,  from  (5.27) 

0       _  1-11     7]^    -      i-oSU  7^ 

tr    -  .  (5.31) 

Besides  the  trivial  solution,  (5  28)  has  the  following  roots: 

rj  =   06033    7    ->?=/6o)  ,  r]  =     -6  02^06.     ^^.^^^ 

Corresponding  to  the  two  positive  real  roots,  there  are  two  values 
of  q  solved  from  (5.31),  namely, 


Note:   This  phenomenon  actually  corresponds  to  Tsien's 
"lower  buckling  load"  in  the  static  problem.   However,  as  we  have 
seen  in  the  static  case,  the  lower  buckling  load  cannot  be  explained 
without  using  the  "energy  criterion,"  but  there  is  no  such  diffi- 
culty presented  in  the  dynamic  case  when  the  trajectories  of  the 
motion  are  considered. 
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r'  =   ^^^9""   •  (5^34) 

Thus,  when  9  =  0.44,  the  parameter  q  can  assume  either  one  of  the 
above  two  values,  which  will  yield  a  common  root  between  equa- 
tions (5.23)  and  (5.25) ^ 

When  q  takes  the  first  possible  value  given  in  (5.33), 
equation  (5=23)  for  9   =  0.44,  assumes  the  following  form: 

which  has  three  real  positive  roots: 

^,  =  oi^lS^     ,     ^^=0  6o35  ,   7^=  ,.7278-2    (5.35) 

Tested  by  condition  (5.26),  we  find  that  (  ^o.  0)  is  the  unstable 
singular  point.   Also,  as  q  =  q    ,  equation  (5.25)  takes  the 
form: 

which  has  three  real  positive  roots: 
7]     =     2l^ZoZ       , 

and  a  double  root  at  ^ =  0.6033o  (5.36) 

Comparing  the  results  in  (5.35)  and  (5.36),  it  is  apparent  that 
"^  =  "^2  "  0.6033  is  the  common  root  of  (5.23)  and  (5.25),   Conse- 
quently, the  trajectory  defined  by  equation  (5.24),  when  0=   0.44, 
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q  =  q    ,  Is  the  separatrix,  because  it  passes  through  the  saddle 
point:   ^  =  0.6033,  ^  =  Oo   For  any  q  >  q    ,  the  motion  of  the 
system  will  be  unstable  and  the  amplitude  of  the  vibration  undergoes 
a  severe  change.   We  shall  call  this  the  "type  I  buckling."^ 

The  second  possibility  that  (5.23)  and  (5.25)  have  a  real 

—   —  (2) 

positive  common  root  is  investigated  by  taking  q  =  q    ,  defined 

in  (5.34).   In  an  analogous  fashion,  we  find  three  critical  points 
by  solving  the  corresponding  equation  (5.23).   Thus, 

rj^  =    0  0  7^3    '  \-o9U36C      '        rj^-    l-Soj   . 

Among   these   three  points,    according   to    (5.26),      11=     ^2  ~  0.84366 
is    the  unstable   saddle  point.      Roots   of    the   corresponding   equation 
(5.25)    are   found   to   be 

■7=  0.15663,  and  a  double  root:   "7=  1601. 

Here  it  is  shown  that  ^=  1  601  is  the  common  root  of  (523)  and 

(5.25).   The  motion  of  the  system  characterized  by  the  parameter 

-   —  (2) 

q  =  q     will  not  become  apparent  without  examining  its  trajectory 

in  the  phase  plane » 

After  substituting  all  values  for  the  parameter  in  (5.24), 
this  equation  takes  the  form: 


1 
Refer  to  Figure  31, 


(5.37) 
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^^  -       r][   (0)^663  -'7)('/.6o|  -  ^f  ]      , 
or  equivalently, 

I   =  lJrj(o)^Q(>-5-ri)(i(,ol-r)f       , 

which  is  the  equation  for  the  trajectory  of  the  motion.   Because 
of  the  branch  points  presented  in  (5.37),  certain  interesting 
behavior  of  the  motion  could  occur  that  may  dominate  the  buckling 
behavior  of  the  shell,  vhich  can  be  obtained  from  the  usual  engi- 
neering tests.   Physically,  '^    is  the  displacement  of  the  motion, 
and  ^  is  the  modified  velocity.   From  equation  (5.37),  it  is  seen 
that  in  the  interval  0.15663  <    "^   <     1.601,  no  real  ^  value  exists; 
i=e.,  no  physical  interpretations  correspond  to  the  motion  of  the 
system  when  '^  falls  into  this  interval.   However,  when  ^=  1.601, 
the  value  of  ^  exists  and  the  system  does  have  meaning.   This  is  to 

say  the  single  point  (  ^=  1  =  601,   ^  =  0)  is  admitted  to  be  one  of 

2 
the  points  on  the  trajectory  of  (5.37),    In  the  real  case,  however, 

the  motion  is  not  defined  in  the  interval  0.1566  <   i^    ^     1.601.   It 

-     -   C) 

has  to  be  explained  that  under  the  action  of  q  =  q     ,  a  sudden 

"iump"  of  the  maximum  deflection  may  occur  as  soon  as  V  reaches 

•'     ^  ■'  max 

the  value  of  0.1566.   The  motion  of  the  system  is  stable  in  the 


This  term  is  borrowed  from  functions  of  a  complex  variable. 
^See  Figure  32. 
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sense  of  the  stability  theorem,  but  buckling  of  the  shell  may  occur 
because  of  the  sudden  and  severe  change  of  deflections;   We  shall 
call  this  situation  the  "type  II  buckling,"  and  the  corresponding 
load  q     as  qcr   ■> 

For  all  other  possible  values  of  d  ,    there  are  also  two 
typas  of  severe  change  of  deflection,  and  the  corresponding  loadings 
which  cause   these  changes-   All  of  them  can  be  found  in  a  similar 
way  as  discussed  in  this  exaraplcc   Results  are  shown  in  Figure  33. 

It  is  obvious  that  the  type  I  buckling  corresponds  to  the 
loss  of  stability  of  the  shell  based  upon  a  classic  criterion  in 
the  static  loading  problems.   Corresponding  to  a  least  applied 
pressure,  a  9   =  0.53  is  found,  and  also  the  corresponding 
critical  loads  are  found  to  be:   q^r    ~  0.665,  q      =0.49. 
The  critical  load  q^^j.    is  seen  to  be  19  per  cent  less  than  the 
static  critical  load  found  in  equation  (5.1A)o 

B,   A  Statistical  Approach  to  the  Determination  of  o^r 

The  motion  of  a  shell  described  by  (5.19)  and  initial  condi- 
tions in  (5  20)  may  represent  the  condition  of  the  shell  under  the 
action  of  a  blast  from  a  -jhock  tube,  with  a  pressure  profile  shown 
in  Figure  6,  Part  B.   Experimental  tests  on  this  problem  show  (20) 


1 
A  linear  problem  of  a  one-degree-of- freedom,  mass-spring 

system  under  such  a  loading  was  formulated  in  (28) « 
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that  the  size  of  the  buckling  region  is  not  unique.    It  depends 
upon  the  shock  velocity,  shell  material,  etc.   We  also  have  reason 
to  believe  that  the  size  of  the  buckling  region  under  a  dynamic  load 
depends  also  upon  all  probable  disturbances  and  imperfections  which 

have  been  shown  to  have  much  more  influence  on  the  shell  buckling 

2 
in  this  case  than  they  do  in  a  static  load  case.     Because  of  the 

non-uniqueness  of  the  size  of  the  buckling  region  and  the  higher 
sensitivity  of  the  shell  buckling  stress  to  all  possible  disturbances, 
we  may  consider  the  size  of  the  buckling  dent  as  a  "random  variable." 
With  the  lack  of  information  from  the  "sample,"  i.e.,  the  experi- 
mental test,  we  shall  set  up  some  plausible  assumptions: 

(i.)   Assume  that  the  number  of  the  facts  which  influence  the 
size  of  the  buckling  region  is  large. 

(ii.)   Assume  that  all  the  facts  (or  disturbances)  can  be  con- 
sidered as  a  set  of  independent  impulses, 

(iii.)   These  independent  impulses  act  on  the  shell  in  order,  i.e., 

3 
one  after  another. 

(iv.)   The  change  of  the  size  of  the  buckling  region  caused  by 

each  impulse  is  proportional  to  the  magnitude  of  the  in^ulse  and 

directly  proportional  to  the  momentary  size  of  the  buckling  region. 


1 

On  the  other  hand,  static  tests  usually  show  a  unique  size  of 

the  buckling  dent. 

2 
For  example,  reference  (12) . 

3 
For   simultaneous  actions,  we  can  simply  consider  them  as 

a  single  impulse. 
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If  the  above  assumptions  are  admitted,  then  the  effect 
on  the  size  of  the  buckling  region  caused  by  various  random 
disturbances  has  the  very  same  nature  as  the  effect  on  the  size  of 
a  certain  organ  caused  by  a  set  of  impulses  in  the  problem  of 
biological  science.   It  has  been  discussed  in  reference  (10)  that  -j 
such  a  random  variable  has  a  log  normal  distribution;   i.e.,  the 
size  of  the  buckling  dent  region  can  be  considered  as  a  random 
variable  having  log  normal  distribution. 

As  it  has  been  defined  in  (5.18),  the  size  of  the  buckling 

2 
dent  is  proportional  to  1/0  ,   and  as  discussed  previously,  6  can 

take  any  value  in  a  certain  interval  defined  in  (5o30).   Therefore, 

we  can  define  a  new  variable  X  , 


O  <  X  <  c^ 


(5.38) 


which  represents  the  size  of  the  buckling  region.   x  has  a  frequency 
function  in  the  form  of 

(5.39) 


Refer  to  Figure  30= 


2 
This  is  assuming  that  R^  is  taken  as  a  constant  value. 
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The  j-th  moment  about  the  origin  (2)    has  the  form: 


(5.40) 


Therefore,  the  expectation  of  x,  i.e.,  the  mean  of  the  distribution 

is 

r  ** 

E[x]  =    X  dA   =  exp  [m   "^  jS^]  , 

^o  (5.41) 

It  is  seen  from  (5.41)  that  the  mean  value  of  the  distri- 

2 
bution  depends  upon  two  parameters:  /-l  and  S  ,  which  are  respectively 

the  mean  and  variance  of  the  distribution  of  the  random  variable  y, 

'^  =    In  X   , 

2 
which  has  a  normal  distribution;  i.e.,  y  is  N(yU,  S  ).   Again,  yU 

2 
and  S   values  can  be  determined  from  the  sample.   Here,  it  is  an 

arbitrary  but  naturel  assumption  that  the  variable  y  =  In  x  has  a 

standard  normal  distribution,  i.e.,  y  has  N(0,  1),  where 


dji  -L 

Therefore,  we  obtained  from  (5.41)  the  value: 


Ij  -  jk  -^[-i>^]  .    -"'»^ 


E[.]   =    e°' 


XwAN   =    e°'=   l-6ir 


(5.42) 
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From  (5.38),  we  can  obtain  the  mean  value  of  Q    by  use  of  (5.42); 
thus 

1 


^m     "^   /o-729 


(5.43) 


If  Ojjj,  the  most  probable  value  of  6  based  on  the  above 
assumptions,  is  taken  to  be  the  Q^^   which  defines  the  size  of  the 
buckling  region,  then  a  similar  analysis  as  in  the  previous  example 
gives  the  corresponding  critical  loads  for  the  type  I  and  type  II 


buckling  as  follows: 

"  cr 

(5.44) 

-  (1) 
From  the  above  equation,  it  is  seen  that  q^^.    is  about 

94  per  cent  of  the  static  critical  pressure  q   given  in  (5=14), 

-   (2) 
while  the  type  II  buckling  pressure  q      is  about  50  per  cent  of 

the  static  lower  buckling  load:  q  =  0.35. 


These  values  may  be  obtained  directly  from  Figure  33. 
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4.   Buckling  under  a  Pure  Impulse 


Following  a  similar  argument  presented  in  Chapter  II  for 
shallov  spherical  shells,  we  hav*.'  the  study  of  the  response  of  a 
complete  spherical  shell  under  a  pure  impulse  in  the  form  of  Part  A 
of  Figure  6,  resulting  in  the  solution  of  the  following  differential 
equations: 

ij   +  [(276  t  J  02^    d^)   7  -  5  3)  f    -(■  2.lo^    Tj^   j 


:=    0     , 

which  is  obtained  from  equation  (5.1y)  by  dropping  the  forcing 
function,  and  together  with  the  initial  conditions: 

77(0)  =  0        ,  ^(0)  =    ^  (^Qf)       , 


(5.45) 


(5.46) 


which  is  obtained  by  considering  the  conservation  of  linear  momentum 
at  the  initial  time;  this  has  been  discussed  previously  in  Chapter  II. 
In  both  of  these  two  equations,  once  more  the  notations  in  equa- 
tion (2  4)  and  (2.5)  for  the  nondimensional  physical  quantities 
have  been  employed.   By  comparing  equation  (5.43)  and  (2,7),  we 
found  that  the  free  oscillation  of  a  complete  spherical  shell  in  the 


Refer  to  Figure  6,  Part  A. 
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region  of  buckling,  or  in  the  boundary  layer  region,  as  one  may  say, 
has  the  same  nature  as  a  shallow  spherical  shell  clamped  along  a 
circular  boundary.   The  only  difference  between  the  two  systems  (2o7) 
and  (5.45)  Is  that  the  size  parameter  Q    is  being  considered  as  a 
variable  in  the  latter  case. 

The  stability  property  of  system  (5.45)  can  be  discussed  in 
an  identical  manner  as  that  performed  for  equation  (2.7).   We  have 
found  that  the  loss  of  stability  of  system  (5.45)  is  characterized 
by  the  condition  that  its  maximum  central  deflection  reaches  the 
critical  magnitude  described  in  equation  (2.19),  i.e.. 


7]^^  =      i-2^9US'  -    I  igi^l  yo  19303  -  e^ 


(5o47) 


The  0  value  which  characterizes  the  buckling  of  a  complete  shell, 

i.e.,   0   ,  may  be  determined  by  two  postulations.   Both  of  these 
cr 

two  postulations  have  been  used  in  the  last  two  examples  without 

particularly  mentioning  that  they  are  based  upon  a  common  concept 

that  the  development  of  the  buckling  region  is  a  continuous  process. 

In  the  first  postulation,  we  have  considered  that  the  size  of  the 

buckling  region  is  a  dependent  parameter  of  the  critical  load,  hence 

the  instability  of  the  shell  is  characterized  by  the  load  and  Q 

J  -'  cr 

is  determined  by  yielding  a  smallest  buckling  load.   In  the  second 
postulation,  on  the  contrary,  we  have  considered  the  development  of 
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the  most  possible  size  of  the  buckling  region  characterized  the 
instability  of  the  shell.   If  a  continuous,  rather  than  random, 
development  of  the  size  of  buckling  dent  is  considered,  then  the 
buckling  of  the  shell  is  detemined  by  forming  a  smallest  size  of 
the  buckling  region  or  forming  the  thinnest  boundary  layer.   It  is 
apparent  if  the  buckling  load  is  a  monotonous  increasing  function 
of  the  buckling  dent  size,  or  vice  versa,  then  the  two  postulations 
vill  yield  the  same  critical  conditions  for  the  shell    It  seems 
that  the  buckling  under  a  static  load  follows  the  first  postulation 
very  well,  while  for  the  buckling  under  a  dynamic  load  the  second 
postulation  will  give  a  better  solution. 

By  consulting  the  result  presented  in  Figure  10,  we  can  see 
that  the  formation  of  the  smallest  buckling  dent,  corresponding  to 
the  largest  Q    value,  is  equivalent  to  yielding  a  least  value  of 
the  critical  impulse.   Therefore,  the  largest  value  of  9    in  equa- 
tion (2.1A)  for  the  limit  of  buckling  of  shallow  shells,  is 
the  critical  6  for  the  loss  of  stability  of  a  complete 

spherical  shell,  i.e., 

(5.48) 

The  corresponding  critical  deflection  is  obtained  from  (5.47). 

^ct     =        '•2^9  .  (5.49) 

The  equation  of  the  phase  plane  trajectories  of  system  (5  45)  will 
be  the  same  as  the  system  in  (2.7),  because  they  have  the  identical 
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fonn„   Therefore,  the  constant  C  which  determines  the  equation  of 
the  separatrix  in  the  phase  plane  of  system  (5.45)  will  have  the 
identical  form  as  equation  (2ol8),  i.e., 

Cj  =    111    r^l^    -     I  o^Jg-  7j^ 

where   "O   takes  the  value  in  (5o49)o   Hence,  we  have 
cr 

Cs   =    0■8gU^^        .  (5.50) 

The  critical  dimensionless  impulse  is  determined  from  the  equation: 

M")  =    f  (».,C)  =  /^     , 

that  is: 


f  = 


(5.51) 


Substituting  (5. 48)  and  (5.50)  into  (5.51),  we  obtain  the  critical 
impulse  for  a  complete  spherical  shell: 

I^,^   =     ).55U5li  .  (5.52) 

As  in  the  last  section,  let  us  consider  9  as  a  random 

variable  and  assume  the  variable: 

I      I 
^      ~  Q    ~       oliil  ,  0  <  0  <  0  iiU 

<^<^<     ^     '  (5.53) 

2     1 

has  a  log  normal  distribution  with  /•  =  0,  S  =1.   Following 


Refer  to  equation  (5.41). 
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the  procedure  performed  in  the  last  example,  from  (5.41)  to  (5.42), 
we  have  the  expectation  for  6   ,  i.e.,  its  mean  of  the  population: 


*^'"        J  6^  +  '/OULt-  —  0  2^U  .  (5.54) 

From  (2.14),  (2.18),  and  (5.51),  we  found 

^,,  =    0.^3   .     &  i;';=  2.oU.  (5.55) 

This  is  about  24  per  cent  higher  than  the  result  in  (5,52), 
obtained  without  considering  0  as  a  random  quantity. 

5.   Conclusions  and  Discussions 

Static  and  dynamic  buckling  of  a  complete  spherical  shell  has 
been  analyzed  based  upon  an  axisymmetrical  theory  and  a  first- 
approximation  governing  equations  obtained  from  a  new  postulation 
of  the  buckling  mechanism. 

A  static  critical  pressure  for  the  inward  snap  was  obtained 
and  has  a  magnitude  about  68  per  cent  of  the  classic  buckling  load 
given  by  Zolley;  while  an  outward  snapping  pressure  was  about  the 
order  of  the  well-known  "lower  buckling  pressure"  by  Tsien,  i.e., 
about  26  per  cent  of  Zolley' s  result.   Both  of  these  two  critical 
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pressures  are  in  fairly  good  agreement  with  recent  experimental 
results  given  by  Thompson  and  Krenzke  (33)  and  (19) . 

A  dynamic  buckling  problem  has  been  studied  under  the  usual 
simplification  introduced  in  shallow  shells;  i.e.,  the  inertia  forces 
in  the  middle  surface  are  neglected.   Based  on  this  simplification, 
a  single-degree-of-f reedom  system  was  obtainedo   Results  obtained 
from  this  analysis  indicate  that  the  "size  of  the  buckling  region" 
seems  to  characterize  the  dyanmic  buckling  of  a  complete  spherical 
shell.   In  contrast,  the  static  buckling  is  characterized  by 
the  "smallest  critical  pressure."   In  accordance  with  the  random 
nature  of  the  size  of  the  buckling  region  under  a  dynamic  loading, 
we  have  assumed  a  log  normal  distribution  for  this  random  variable. 
A  statistical  approach  based  on  an  assumed  probability  distribution 
gave  a  critical  pressure,  for  the  shell  under  a  dynamic  load  in 
the  form  of  a  step  function,  with  magnitude  about  94  per  cent  of  the 
critical  static  pressure  (based  upon  (^P  static  "  0.821).   It  has 
been  pointed  out  in  the  analysis  that  a  second  type  of  buckling  may 
occur  due  to  a  "jump"  of  the  maximum  central  deflectiono   Based  on 


Note:  The  critical  pressure  predicted  by  this  theory  is 
very  close  to  the  result  of  Krenzke' s  experiment  but  somewhat  lower 
than  the  result  given  by  Thompsono  We  think  that  perhsps  this  is 
mainly  due  to  the  difference  in  Poisson's  ratios.  The  model  used 
by  Thompson  was  polyvinyl,  which  has  a  Poisson's  ratio  i>  =  0.48, 
while  Krenzke' s  aluminum  model  has  ^=  0.3,  which  is  used  in  this 
analysis. 
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the  statistical  approach,  we  have  obtained,  corresponding  to  the 
outward  snapping,  the  pressure  with  magnitude  about  50  per  cent  of 
Tsien's  "lower  static  buckling  pressure."   If  we  define  a  quantity: 

(  (fci   static 
(Dynamic  buckling  factor)  =  , 


^  Ucr  /  dynamic 

(5.56) 

then,  according  to  (5.A4),  a  DoB.F.  =  2.14  is  obtained.   A  recent 
experimental  test  (20)  reported  an  average  value  of  D.B.Fo  =  2.18, 
which  <ias  obtained  by  testing  20-inch-diaraeter  aluminum  spherical 
shells. 

Buckling  of  a  complete  spherical  shell  under  a  pure  impulse 
has  also  been  studied.   The  critical  condition  was  determined  in 
such  a  manner  that  the  loss  of  stability  is  characterized  by  forming 
the  smallest  region  of  buckling.   In  this  way,  we  obtained  a  non- 
dimensional  critical  impulse  for  a  complete  spherical  shell. 

The  legitimacy  of  neglecting  the  inertia  forces  in  the  middle  surface 
was  not  justified;  it  remains  an  open  question.   It  seems  to  the 
author  that  this  might  introduce  a  larger  error  than  in  the  case  of 


Note:   This  value  is  obtained  by  taking  (icr^ static  ~  0«32 
given  in  equation  (5.15). 
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shallow  shells,  because  the  membrane  stresses  play  important  roles 
in  the  earlier  stage  of  deformation  of  a  complete  shell.   However, 
strange  as  it  may  seem,  we  found  in  most  dynamic  considerations  of 
this  problem,  such  as  the  one  in  reference  (A),  these  inertia  forces 
were  also  neglected,  even  when  a  membrane  shell  theory  was  employed. 

In  a  paper  on  "Dynamic  Criteria  of  Buckling"  (13) ,  in  1949, 
Hoff  proposed  a  test  of  stability.   He  suggested  that  stability 
investigations  be  carried  out  by  assuming  a  probable  disturbance 
of  the  equilibrium  and  calculating  the  ensuing  displacements  from 
the  dynamic  equations  of  motion.   The  nature  and  the  magnitude  of 
the  disturbance  must  be  established  from  a  statistical  investigation 
of  the  conditions  under  which  the  structural  element  or  the  part  of 
machinery  will  be  used.   The  safety  of  the  system  can  be  safeguarded 
if  it  is  made  stable  for  all  disturbances  which  have  a  probability 
greater  than  a  required  minimum.   This  idea,  which  has  not  received 
as  great  attention  in  the  Western  World  as  it  should  have,  has  been 
extended  systematically  by  several  Russian  authors,  especially 
Vo  V.  Bolotin  (5)  »   We  attempted,  in  this  research,  to  establish 
a  suggestion  that  Hoff's  idea  may  have  much  significance  in  the 
dynamic  buckling  problems.   Instead  of  considering  all  probable 
disturbances  that  would  make  the  nonlinear  equation  more  complex  and 
be  a  branch  of  research  of  its  own,  we  have  allowed  a  probable  effect 
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caused  by  the  probable  disturbances  on  the  deformation  of  the 
structure,  which  can  be  expressed  in  the  form  of  a  parameter.   The 
most  probable  critical  condition  of  the  structure  is  determined 
by  a  specific  value  of  this  parameter  obtained  in  a  probability 
sense.   Experimental  evidence  shows  that  the  shell  buckling  under 
a  dynamic  load  (rapid)  is  different  in  character  from  those  under 
a  static  loading  by  forming  non  unique  buckling  patterns,  e.g., 
in  spherical  shells  the  size  of  the  buckling  dent  varies;  in 
cylindrical  shells  the  number  of  dent  varies  (37) .   Of  course, 
these  quantities  seem  dependent  more  or  less  upon  other  physical 
quantities  such  as  pressure,  speed,  materials,  etc.;  they  are  also 
influenced  by  external  disturbance  and  imperfections.   Studies  in 
(1_)  ,  (12)  ,  and  (18)  show  that  the  shells  are  much  more  sensitive 
to  the  disturbances  and  imperfections  under  a  rapidly  applied 
dynamic  load  than  under  a  load  applied  gradually  as  in  a  static 
sense.   It  seems  more  desirable,  as  shown  in  this  investigation, 
that  the  effects  due  to  disturbances  on  the  buckling  of  the  shell 
be  assumed  in  a  random  nature.   Therefore,  from  the  designer's 
point  of  view,  a  further  study  based  upon  the  theory  of  Hoff  and 
Bolotin  would  be  very  interesting,  as  far  as  obtaining  a  safe  load 
is  concerned. 


CHAPTER  VI 

A  CONCLUSIVE  REMARK 

A  dynamic  buckling  criterion  which  is  independent  of  the  type 
of  loading  and  geometrical  shape  of  the  shell  structure  has  been 
proposed.   This  criterion  is  founded  on  a  comparative  basiSo   The 
critical  condition  of  the  structure  is  determined  from  its  dynamic 
characteristic .load-characteristic  deformation  relation,  in  compari- 
son with  the  critical  conditions  of  the  static  loading  and  free 
oscillation  of  the  same  structure  for  conservative  and  non- autonomous 
systems.   Examples  in  Chapter  II  shows  that  the  critical  condition 
determined  by  this  criterion  is  identical  with  that  determined  by 

the  stability  theorem  of  nonlinear  mechanics   and  agrees  well  with 

2 
the  results  obtained  by  other  means. 

A  unified  study  has  been  made  without  particular  mention. 

As  we  know,  the  general  property  of  the  stability  of  motion  of  a 

dynamic  system  can  be  directly  determined  by  finding  the  existence 

of  a  certain  positive  definite  function--the  Liapunov  function.   In 

case  that  the  system  is  autonomous  and  conservative,  the  Hamitonian 


This  name  was  used  by  Davis.   Refer  to  (9),   p.  317. 
^Refer  to  Figures  10  and  14. 
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is  known  to  be  the  Liapunov  function  of  the  system.   Therefore,  the 
value  of  the  parameter  determined  by  the  critical  condition  of  the 
Liapunov  function,  or  in  this  case,  the  Hamitonian,  is  by  no  means 
a  less  accurate  one  than  that  which  would  be  obtained  from  the  study 
of  dynamic  response;  further,  it  is  much  more  formal  than  by  the 
latter  method.   However,  in  some  other  cases,  when  the  system  is 
nonautonomous,  there  is  no  general  method  to  find  the  Liapunov 
function,  hence  the  technique  used  to  determine  the  critical  parameter 
equivalent  to  the  energy  method  might  not  be  possible  for  such  caseSo. 
Furthermore,  where  the  dynamic  system  is  known  to  be  unstable,  then 
the  energy  method  will  not  lend  itself  to  the  determination  of  the 
critical  condition  of  the  structureo   It  is  seen  from  the  examples 
that  by  a  phase-plane  topological  method  the  critical  parameters  can 
be  obtained  in  a  very  simple  manner,  which  is  equivalent  to  that  by 
an  energy  method;   however,  it  gives,  in  addition,  a  general  quali- 
tative picture  of  the  motion  of  the  system.   Furthermore,  established 

2 
methods  are  available  to  construct  the  time  history  from  the  trajec- 
tory of  the  system,  hence,  the  response  curves  can  also  be  obtained 
in  a  relatively  easier  way  if  it  is  desired. 


Refer  to  reference  (15) 

2 
Refer  to  reference  (42) 
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A  general  relation  between  the  critical  load  and  geometrical 
parameter  in  the  case  of  a  suddenly  applied  constant  pressure  was 
found.   It  is  seen  from  Figure  12  that  this  relation  has  approx- 
imately the  same  nature  as  its  static  counterpart.   The  buckling 
pressure  is  about  75  per  cent  of  the  critical  static  pressure,  and 
the  deflection  of  the  shell  at  the  time  of  buckling  is  seen  to  be 
about  1.6  times  the  static  critical  deflection. 

Transitional  instability  of  the  shallow  shell  under  the 
action  of  a  uniform  pressure  linearly  increased  with  time  was  also 
studied.   Results  are  obtained  by  the  use  of  the  proposed  buckling 
criterion  and  are  summarized  in  the  corresponding  curves.   The  region 
of  the  transient  oscillation  is  stretched  along  with  the  increment 
of  the  rate  of  the  applied  load   This  same  feature  was  also  found 
for  a  cylindrical  shell  by  Vol'mir  and  Agamirov  ►   The  initial 
imperfections  were  not  considered  in  this  study.   They  are  expected 
to  have  significant  influence  when  the  rate  of  load  is  very  high. 

It  is,  therefore,  suggested  that  a  further  research  based 
upon  both  higher  axisymmetric  and  nonsymnetric  modes  of  deformations 
and  also  initial  imperfections  be  studied.   The  effect  of  these 
quantities  to  each  loading  condition  may  then  be  decided  after 
comparison  with  the  present  theory. 

Transitional  dynamic  instability  of  a  complete  sphere  was 
studied  by  neglecting  the  longi tudingly  inertias.   This  simplifica- 
tion was  not  justified.   However,  according  to  Bolotin  (4^),  the 
errors  caused  by  introducing  this  simplification  may  be  expected 
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to  be  small,  if  the  frequency  of  the  external  pressure  is  siaall 
as  compared  with  the  natural  frequency  of  the  longitudinal  vibration 
of  the  system.   The  problems  which  have  been  considered  in  this 
paper  were  certainly  the  cases  which  the  above  argument  validates. 
It  seems  to  the  author  that  Bolotin  arrived  at  the  conclusion  in 
(4)  that  in  the  study  of  the  nonlinear  oscillation  of  a  thin 
clastic  sphere,  one  has  to  consider  the  shell  as  a  whole.   However, 
as  far  as  the  buckling  of  the  shell  is  concerned,  we  feel  that  it  is 
satisfactory  to  consider  just  a  segment  enveloped  by  the  boundary 
of  the  buckling  region.   Yet,  when  the  shell  is  under  the  action  of 
a  periodic  external  force,  we  may  havr  to  consider  the  whole 
spherical  shell,  because  the  buckling  of  the  structure,  most  likely, 
does  not  occur  in  the  transient  state.   As  well  as  in  the  static 
case,  many  more  experimental  tests  on  the  instability  of  a  complete 
spherical  shell  are  highly  desirable. 
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APPENDIX  I 

THE  STABILITY  THEOREM  OF  NONLINEAR  MECHANICS^ 

The  Stability  Theorem 

In  the  neighborhood  of  a  singular  point  (  '^q     |.  )  of  the 
following  system  (A.l),  the  stability  characteristics  of  the  system 
are  determined  by  the  roots  of  its  characteristic  equations. 

The  Characteristic  Equation  of  the  System 

Given: 

i    =      Pc; ,  M 

^    =     ^"'-^^  (A.l) 

vhere  P  and  Q  are  analytical  functions  of  ^7  ^^'^   ^  >    therefore, 
they  may  be  expanded  in  the  form  of  a  Taylor  series  in  the  neigh- 
borhood of  (  ''J       4o)'   Let  us  make  the  following  transformation: 

^        =  ^  +  >.  (A. 2) 

Therefore,  the  study  of  the  original  system  in  the  neighborhood 


The  major  part  of  this  section  is  taken  from  reference  (9), 
pp.  317-318„ 
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of  the  point  (  ^   ,     ^  )  is  equivalent  to  the  study  of  the  new  system 
in  the  neighborhood  of  the  origin  (x  =  y  =  0) „   By  the  transforma- 
tion (A. 2),  (A.l)  can  be  written  as 

In  the  case  that  (  ''7q,  ^  )  is  the  singular  point  of  the  system  (A.l), 
i.e.,  P(  r^Q,     ^q)  =  Q(  i7o'  ^o^  "^  ^'  ^^  ^^^^  ^^^   above  system 
in  the  following  form: 


(A. 3) 


This  is  a  nonlinear  system  equivalent  to  (A.l),  where  f , (x  y)  and 
f2(x,  y)  are  nonlinear  functions  in  x  and  y  of  higher  orders  than 
the  second.   As  we  have  mentioned,  the  study  of  the  motion  in  the 
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neighborhood  of  the  singular  point  (  ^  ,  H,  )  of  system  (A.l)  is 

o    o 

equivalent  to  a  study  of  the  motion  in  the  neighborhood  of  the 
origin  (0,  0)  of  the  system  (A.3)o 

We  shall  call  the  following  equations  the  linear  approx- 
imation of  the  system  (A. 3),  thus. 

The  characteristic  equation  of  (A. 4)  reads: 


(Ao4) 


A  (A)  = 


=     0.     (A. 5) 


We  have  the  following  theorem  regarded  as  the  application  of  (A. 5) 

to  the  nonlinear  system  (Ao3): 

Theorem.   If  the  coefficient  matrix  of  (A. 4)  is  non-singular;  i.e.. 


^. 


0, 


(A.  6) 


and  both  f,  and  f  are  continuous  and  satisfy  the  following 
nonlinearity  condition,  thus 


^Refer  to  reference  (_29) ,  pp.  175-178  and  p,  132. 
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1  ^,^^'V\   +  li;(=t'V\ 

Lira  . 

(A.  7) 


Then  the  stability  character  of  the  critical  point  at  the  origin 
(0,0)  of  the  nonlinear  system  (A. 3)  is  the  same  as  its  linear 
approximation  (Ao4)o   The  stability  of  (A. A)  is  determined  by  the 
characteristic  equation  (A.5)o 

As  an  example,  let  us  consider  (A.l)  as  the  system  (2.8). 
Subjected   to  the  transformation  (A. 2),  (Ao3)  is  in  the  form: 

i  =   0  4   ^   , 

^  =   -[C2  76  +  5  025  eM  -  10  U  r)^   +  C32U  rj^  J    ^ 

It  is  seen  from  the  above  equation  that 
j-,   i^,*^)       =         0    , 

Then  the  nonlinearity  condition  (A. 7)  is: 

Lim  1  (5-3)  t  6  32U  V^)  -   2io2   X  /  I  X^j 

[HI,  l^]->o  ,     ,  .  ^ 

hence,  the  nonlinearity  condition  is  satisfiedo 
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Therefore,  the  characteristic  equation  (A. 5)  can  be  used  to  study 
the  nonlinear  system  (A, 4),  and  it  has  the  form: 


:=   0 


-X  ; 

-[(276  ■^  302^   eM  -  10  62  ^^  +  (,-32k   ^/ ]       "A 

which  is  seen  to  be  equation  (2„12). 

The  stability  characteristic  is  determined  by  the  following 
rules: 

A.  When  both  roots  of  the  characteristic  equation  (A. 5) 
are  purely  imaginary  values,  the  singular  point  is  a  "center,"  and 
the  motion  in  the  neighborhood  of  this  point  is  stable. 

B.  When  both  roots  of  (A. 5)  are  real  but  differ  in  sign, 
the  singular  point  is  a  "saddle  point"  and  the  motion  in  the  neigh- 
borhood is  unstable,  etc. 
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THE  EQUATION  OF  MOTION  FOR  TRANSVERSE  VIBRATION  OF 
SHALLOW  SPHERICAL  SHELLS 


It  is  well  known  that  for  a  thin  shallow  elastic  shull, 
the  following  important  simplifications  are  permitted  (21) : 
Ao   The  effect  of  transverse  shear  is  neglected, 
B.   Relations  between  moments  and  changes  of  curvature 
are  kept  the  same  as  in  the  theory  of  plates.   Therefore,  the 
strain  energies  due  to  membrane  and  bending  deformations  are  the 
same  as  in  the  case  of  a  plate,  thus 


2(i-V^) 


J  V 


(A.  8) 


JJJV   "• 


(A.  9) 


where  a.  and  b.  are  the  first  and  second  invariants  of  the  displace- 
ment;  their  definitions  will  be  given  later. 

Referring  to  Figure  1,  in  polar  coordinates,  we  have  the 
strain-displacement  relations: 


Refer  also  to  reference  (26) 


124 


125 


(A. 10) 

^ee  =  -[j^r   +  72  Wee]  n 
■^re  =  -  [  7  "^re  -  72  ^0  ]  ^ 

(A. 11) 

In  (A. 10)  and  (A. 11),  r,  Q     are  the  conventional  polar 
coordinates;  ''1  is  the  variable  in  the  direction  normal  to  the 
middle  surface  and  measured  along  the  thickness  of  the  shell; 
u,  V,  and  w  are  the  displacements  in  the  directions  of  r,  0  ,  "t 
respectively;  z  is  the  equation  of  the  surface  of  the  shell,  i.e., 
the  vertical  distance  measured  from  plane  TT  to  the  surface  Z  as 
indicated  in  Figure  1.   It  is  conventional  to  take  the  surface 
of  a  shallow  spherical  shell  as  a  paraboloid,  which  has  the 
following  expression 

Z  s:   Zo  (  ^/To)^  (A.  12) 

^r    --  ^o  (A. 13) 
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For  a  laterally  loaded  shallow  shell,  let  us  consider  that  the 
nonlinear  terms  in  u.and  v  and  their  derivative  with  respect  to 
r  and  9  ere  small  compared  with  the  transverse  displacement  of  the 
structure.   Therefore,  by  introducing  (A. 12),  (A. 13),  and  the  above 
simplifications,  (AdO)  reads: 

606=   ?(u+  Ve)  +  i  (7^V0)  . 

(A. 14) 

By  using  (A. 14)  and  (A. 11),  we  have  the  invariants  a^^  s  and 

t 

b^  s  in  (A. 8)  and  (A. 9)  in  the  following  form: 

+  vWr IVg  -  -Uq  wr  kVe  )  -  J  (  v' +  ^e  )]  I 
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b,     = 


b,     = 


2-     s  2        _  12-1 


Hence,    (Ao8)    and    (A. 9)    turn  out   to  be  of   the   following  form: 


(A. 15) 


9,= 


ad-i**)    J. 


Iv 


(  a;   +1  W^     +  tifl^r)  -    -^   (^ 


+  l/g   )  iV^     +    f  (2i(U^   +  2CVVe    -t-  t(w/  +  V^  VJr^) 


-   7h.    ^r^9 


'    .2  .  X 


+  72  (m^  ti;/  +  2uVe  ^i^^^^) 


-  xig  vv^  iV0   -   £  v^  -  J  ^e   )  J  I    r  dr  de  dn. 

(A016) 
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E  fff         f         2  2v'  1^ 


kVr.Wg 


+  ^^      4      -,    iV.Wge     +     ^, 


w, 


90 


+  (i-W     -  ^'^'•9      +     U_vJe_^     _    2^ 


2    ^ 


r^    J 


•n^  ro/r  dB  dn  . 


The  work  done  by  the  external  surface  forces: 


(A. 17) 


W  = 


^  w  r  dY  dd 


(A. 18) 


The  kinetic  energy. 


T    = 


2  f  P  ^tt  1  1^  c(©cjr  dn  , 


(Aa9) 


where  p  is  the  mass  density  of  the  shell  material.   In  (A. 19),  we 
have  also  neglected  the  longitudinal  inertias. 
Let  us  designate 

F=     a,  +   92+W-T  (A. 20) 

then,  the  Euler-Lagrage' s  equations  have  the  following  form: 


da      atUu/'     dtiau;./      dev^ug/ 


=  0 
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if-   -(^]    -  i/i^i    -    ^(^\ 

-       0 

> 


n2 


(  c^WflflJ 


^e'  ^  c9Wg3 


=       0 


(A. 21) 

After  substituting  equations  (A. 16)  to  (A. 19)  into  (A. 20), 
and  then  to  (A. 21),  we  have  three  equations  of  motion  in  terms  of 
the  displacements  u,  v,  and  w.   However,  it  is  well  known  that  these 
three  equations  can  be  reduced  to  two  equations  by  introducing  a 
stress  function  (p  ,  which  has  the  definition  in  the  following 
equations: 


0^. 


n    =         ji      +       T^<^6 
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(A. 22) 


It  was  found  by  using  this  function  that  the  longitudinal 
displacements  which  satisfy  the  following  compatibility  condi- 
tion will  satisfy  the  two  equations  of  motion.    The  compatibility 
equation  has  the  form: 

The  third  equation  of  motion  in  the  transverse  direction, 
in  case  the  external  pressure  is  uniformly  distributed  over  the 
shell  surface,  was  found  to  be: 

^  ^  (A. 24) 


Actually,  these  are  the  equilibrium  equations,  because  the 
inertia  terms  were  neglected. 
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Equations  (A. 23)  and  (A. 24)  have  the  same  form  as  found  by  Grigolyuk 
in  the  paper  given  in  (52) ,  where  a  direct  transfonnation  from 
Marguere's  equation  (21)  was  performed.   In  both  equations  (A. 23)  and 
(A. 24),  the  operator  ^   has  the  following  form,  thus 


(A. 25) 


Eh-^ 
and,  D  =   .  „  ,. 7-r—  is  the  bending  rigidity  of  a  shell  elemento 

In  the  axisymmetric  theory,  we  consider  the  displacement 

and  stresses  are  functions  of  a  single  space  variable  r.   Therefore, 

equations  (A. 23)  to  (A.2j)  assume  the  following  simplified  forms: 


(A. 26) 


D  VV    +  k\-^     it^T^r)    -      ^^   -    4^^ 


R.      .  . 


<J   +   lipW,^      :=       0  , 


(A. 27) 


and 


^r)    -       [  r  ^r(^arjj 


(A. 28) 


Further  simplifications  are  possible  by  using  (A. 28)  and 
integrating  equations  (A. 26)  and  (A. 27)  from  0  to  r;  thus,  we  have 
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(A. 29) 


D 
d 


0   , 


(A. 30) 


where, 


These  are  the  governing  equations  used  in  this  dissertation. 

Let  us  assume  a  first-degree  approximation  mode  for  the 
axisyrametrical  deflection  as  follows: 


IV 


=      <a)f  I-  (VrJ]'     ,  (^.31) 


which  satisfies  the  boundary  condition  of  a  clamped  edge.   With 
equation  (Ao31)  into  (A. 29)  and  subjected  to  the  following 
boundary  conditions  for  the  stress  function  (^  ; 

at  r  =  0,      ^   =   0  , 

dr 

U    =     0     ,         or         {<^rr   -   1  ^^r)     =     0  , 


equation  (A. 29)  can  be  integrated;  it  has  the  following  form: 


(Ao32) 
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6  To 


i^^^l-^ib'^^f-ii) 


Ho     I    (>->')  ^^o'  ^"rj  \Yj     _ 


6 


(A. 33) 


Galerkin's  method  permits  that  equation  (A. 31)  be  an  approx- 
imate solution  of  equation  (Ao29)  and  (A. 30),  provided  the  following 
condition  is  satisfied,  i.e., 

t2 


G 


I  - 


d^     =      0 


(A. 34) 


where  A  is  the  total  circular  area  with  radius  r  ,  and 

o' 


G   = 


-  6['  -f  W'  -  ( 


^J 


6 


-  K'   +  i 


J 


^f,^'f(^j''-^' 


.  */'  -;'] 
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5R.rr,    ^    I-      i-v  J 


2D  2D       ^tt  1    ^ 


2[      +    [      I      , 


where   1   =  — 
r„ 


(A. 35) 


After  performing  the  Galerkin  integral  (A. 34),  we  have  the 
dynamic  equation  of  the  central  deflection  4  ^^  follows: 


12 


i\       ^        ~     P^Jo'l 


E       f^n-znv 


6279 
270U 


30 


-   Ti'     ^ 


prU    6()-2;) 


li  ?■ 


i6 


r^ 


(A. 36) 
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When   i)  =   0.3,    (A. 36)    is    read: 


1^ 


'0    62  _2         ^         Sk3         rll 


Ro  ^o 


11 1 

)6     pf, 


(A„37) 
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APPENDIX  III 
FIGURES 
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Figure  1.   Geometry  and  deformations  of  the  shell 


137 


Figure  2.   Phase  plane  trajectories  and  the  variation 
of  potential  energy 
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T=z. 


Figure   3c      Phase   space   trajectory  and  its  projection 
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Case  I 


Case  II 


Case  III 
Figure  4.   The  dynamic  buckling  criterions 
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Figure  5o   Conmarison  of  axi symmetrical  theories  on  the  static 
buckling  of  shall ov  spherical  shells 
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Figure  60      Dynamic   loadings 
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Buckling  criterion, 
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Figure  9.      Critical   in^ulse  determined  by  buckling  criteriona 
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Figure  10c   Comparison  of  analytical  theories  on  dynamic  buckling 
of  shallow  shells  under  pure  impulse 
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Figure  11 o   The  threshold  of  instability  of  a  shallow  shell 

(with  e  =  0.26)  under  a  suddenly  applied  uniformly 
distributed  dynamic  loading 
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Dynamic  loadings 
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reference  (30) 
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Relation  between  the  critical  pressure  and  the  geometrical 
parameter  9  for  shallov  shells  under  the  action  of 
uniformly  distributed  static  and  dynamic  loadings 
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Figure  13o   Relation  between  the  critical  deflection  and  the 
geometrical  parameter  Q     for  shallow  shells  under 
various  dynamic  and  static  loadings 
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Figure  17.   Upper  and  lower  values  of     Z       as  a  function  of  0 
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Figure  19.   Critical  D.O.L.F.  --  7   vs.  critical  central 
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Figure  20.   Upper  and  lower  critical  deflections  vs. 
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Figure  22^  i       vs.  "(J/  curves  for  various  values  of  Q 
cr 


157 


1 


-  5 


10 


2oO 


_-5 


Figure  23 o   The  projection  of  the  trajectories 

on  the  ^  '   - '^     plane,  when  9  =  0.26, 
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Figure  25.   Buckling  region  of  a  complete  spherical  shell 
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Figure  26o   Load  deflection  curve  for  a  complete 
spherical  shell 
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Figure  27.   Comparison  of  theoretical  result  and 
experimental  tests 
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when  q  •»  q 


when  q  >  q 


Figure  28.   Sketched  phase  plane  trajectories  for  the 
motion  of  a  shell 
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Figure   29.      Graphical   solution  of  equation   (5.29) 


Figure  30 „      Frequency  curves  of  the  normal  and  log-normal 
distributions 
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Figure  32o  'Phase  plane  trajectories  when  9  =  0o44,  q  =  0o35 
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Figure  33.      Buckling  pressure   vs,   buckling   region  parameter   9 
for  a  complete   spherical   shell 
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APPENDIX  IV 

NOTATIONS 

Flexural  rigidity  of  a  shell  defined  by  D  =  Eh^/12(l  -  il^) 

E       Young's  Modulus 

h       Thickness  of  a  shell 

I*      The  dimensionless  inpulse  defined  as:   I*  =  q(^iT) 

N       Distribution  function  of  a  normal  variable 

p       The  critical  pressure  of  a  conplete  sphere  given  by 

2 


the  classic  linear  theory,  p^   = 


V3(l  -    V  h 


q       External  static  or  dynaraic  load  distributed  uniformly 

on  the  shell  surface 
q       Nondimensional  fom  of  the  external  pressure  q  defined 

in  equation  (2.4) 
Q       The  pressure  increasing  rate  defined  in  equation  (4.4) 
Q       The  nondimensional  form  of  Q 
r       Radius  of  a  circular  region  of  the  buckling  dent  of  a 

complete  shell  or  the  radius  of  the  base  of  a  shallow 

shell 

R       Radius  of  the  spherical  shell 

2 
s       Variance  of  a  normal  distribution 

t       Time  variable 
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V  Volume  parameter 

V  Speed  of  propagation  of  the  elastic  waves  in  the  shell 

material,  defined  in  equation  (4.2) 
w(t)     Deflection  of  the  shell  normal  to  its  middle  surface 
X       A  random  variable  defined  in  equation  (5.53) 
y       A  random  variable  defined  as  y  =  In  x 
^       Central  deflection  of  the  shell  or  the  amplitude  of 

w(t)  defined  in  equation  (A. 31) 
z       The  total  rise  of  the  buckling  region  as  shown  in  Figure  1 
7]  The  nondimensional  form  of  C,    defined  in  equation  (2.4) 

Q  A  parameter  to  measure  the  size  of  the  buckling  region 

of  a  complete  shell  or  the  geometrical  parameter  of 

a  shallow  shell  as  defined  in  equation  (2.4) 
/[  Distribution  function  of  a  log-normal  variable 

/*      Mean  of  a  normal  distribution 
1^      Poisson's  ratio 

|,      Nondimensional  velocity       ^    =  '  /  dZ 

P  Mass  density  of  an  element  of  the  shell  material 

(JJl      Radial  stress  in  the  shell  surface 

"J      Dimensionless  time  variable  as  defined  in  equation  (2.4) 
■J*     Dynamic  overload  factor  defined  in  equation  (4.5) 
"vj/      A  nondimensional  quantity  defined  in  equation  (4.7) 
3      Energy  parameter 
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